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16 Abstract 

The subject of this investigation is the propagation of high intensity sound waves through an air-filled porous 
material. The material is assumed (1) to be rigid, incompressible, and homogeneous, and (2) to be adequately 
described by two properties; resistivity r and porosity J2. The resistivity was measured as a function of velocity 
for static flows and found to follow the empirical relation r=ri +r 2 ,u sgn(u). This relation is assumed to apply fo.* 
acoustic signaL as well. Ordinary hydrodynamic .nonlinearity (which leads to slKKk formation) is neglected 
because of the very high attenuatiai in the porous material. The resulting wave equation is still nonlinear, 
however, because of the u sgn(u) term in the resistivity. The equation is solved in the frequency domain as an 
infinite set of coupled inhomogeneous Helmholtz equations, one for each harmonic. An approximate first integral 
formulation of these equations gives relations for progressive waves. The source wave considered is a slightly 
distorted intense tone, that is, a finite- amplitude fundamental accompanied by weak higher harmonics. An 
approximate but analytical solution leads to predictions of excess attenuation, saturation, and phase speed 
reduction for the fundamental component. A more general numerical solution -s o.ed to calculate the propagation 
curves for the higher harmonics. The u sgn(u) nonlinearity produces a cubic distortion pattern; when the input 
signal is a pure tone, only odd harmonic distortion products are generated. Quantitative experiments were 
performed on batted Kevlar*" 29 having porosities in the range f2 =0.94 to 0.98. Measurements were taken in the 
400 to 6200 Hz frequency range for small-signal waves, 900-1500 Hz for finite-amplitude waves. Source levels 
for the finite-amplitude waves extended to 160 dB. Qualitative confirmation of predictions has been obtained in 
all cases, and quantitative confirmation in most cases. The theory is then applied to the problem of propagation 
and attenuation in lined ducts. A vector version of the wave equation is derived along with the corresponding 
inhomogeneous Helmholtz equations. These are used as the basis of a perturbation solution for reflection from a 
porous half-space and propagation in lined ducts. The solution becomes extremely complicated, however, when 
used for high intensity waves in lined ducts. An ad hoc model is therefore presented for the purpose of illustrating 
the gross nonlinear effects on absorption. 
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CHAPTER 1 


INTRODUCTION 


Porous materials have proved useful as sound absorbers in higi- invensity 

environments, such as jet engine inlei ducts, where sound levels routinely exceed 

140 oB. At this level most porous materials exhibit amplitude dependent properties, 

1 2 

i.e., they act nonlinearly. ’ Although it is known that the nonlinearity significantly 
alters the efficiency of the absorber, the physical processes which cause the 
nonlinearity have not been well understood. Moreover, the precise effects of 
nonlinearity are not well documented. There is thus a need to understand more fully 
the high amplitude behavior of bulk porous materials. 

Propagation of sound in porous materials is fundanientally different from 
that in an ordinary fluid. The acoustic flow encounters frictional resistance at each 
.he many gas/fiber interfaces. As a result the sound wave is attenuated, the 
phase speed reduced, and tl.e impedance in<~reased. The severity of these effects 
depends on frequency. For example, low frequency signals pass through tne material 
more by diffusion than by propagation; attenuation and phase speed are both 
proportional toVf . At high frequency, on tli? other hand, the attenuation and phase 
speed reach constant values. Furthermore, because the frictional resistance 
increases with flow velocity, the acoustical properties also depend on amplitude. 
The combination of nonlinearity and diffusion causes some unusual amd interesting 
behavior; it also makes the problem difficult to solve. 

The present study covers high intensity sound propagation in air-filled 
fibrous bulk porous materials. The investigation is both theoretical and experi- 


mental. 


The theoretical work has two hallmarks; (1) the theoretical mode) is very 
simple, and (2) the analysis is done in the frequency domain. The result is a 
relatively uncomplicated theory amd predictions that compare favorably with 
experimental data over a wide range of intensities, frequencies, and input signals. 

Since the number of investigations and treatments of porous material 
problems is vast, it is useful to show the relationship of the present work to that 
which has gone before. Because a tho.'ough review of the literature has already 
been given by Kuntz, ^ we offer aui abbreviated review to establish a succession of 
reference points. 

A. Small-Signal Theory 

3 

Rayleigh was apparently the first to correlate sound absorption in 
porous materials with viscous and thermal losses occurring across the substantial 
internal surface area of the material. He characterized the material as a solid 
riddled with a multitude of narrow tubular channels. Rayleigh combined this model 

4 

with the work of Kii .hhoff who had investigated the effect of viscosity and heat 
conduction on sound traveling in a narrow tube, to create tne first modern theory of 
sound propagation in porous materials. 

Many investigations which followed were based on Rayleigh's pioneering 
5 6 7 8 9 

work. • • ' > The most prominent and thorough work along this line was 
performed by Zwikker and Kosten. They took into account fiber motion, and also 
introduced the concept of the structure factor. 

Lambert ^ ^ discovered that above a certain frequency the gas motion and 
frame motion are indeed decoupled. The decoupling frequency depends on the 
material. For Kevlar*^29, which figures heavily in the experimental phase of this 
work, the decoupling frequency is below 85 Hz for porosities above ft =0.9^. Support 


2 



for Lamberts hypothesis is offered in Chapter 2. H's work increased confidence in 

subsequent theories in which frame motion was neglected altogether. Kuhl and 
1 y 

Meyer were the first to introduce the following two crucial assumptions; (1) the 
dc flow resistance can be u^ed directly in the momentum equation for acoustic 
signals, and (2) frame motion is negligible. Their theory provides the foundation for 
the present study. 

Hersh and Walker elaborated on Kuhl and Meyer's theory. They used 
hydrodynamic flow theory to formulate empirical relations between viscosity, 
porosity, fiber diameter, and resistivity. Heat transfer effects were also accounted 


Other influential linear theories (see, for example, Biot^^ and Lambert*^) 
are comprehensive and generally quite useful. We are not able, however, to derive 
any benefit from them because their inherent complexity precludes extension for 
nonlinear behavior. 


B. High- Amplitude Behavior 

2 

Zorumski and Parrot combined theoretical work with experimental 
observation in their study of +he amplitude dependent impedance of thin porous 
sheets. Most important, they showed that the ac resistivity is not only independent 
of frequency, but that it closely approximates the dc flow resistivity over a wide 
range of particle velocities as well. The equivalence of the ac and dc resistivity 
apparently proves the first assumption of Kuhl and Meyer. 

Kuntz's investigation^ builds on the best aspects of the contributions of 
Kuhl and Meyer and Zorumski and Pa-rot. However, a single error in Kuntz's 
expression for the nonlinear resistivity caused his theoretical predictions to differ 
markedly from his experimental observations. 


3 


There have been other investigations of the nonlinear properties of porous 
materials, but most are based on empirical observation of nonlinearly induced 
phenomena, with no clue about the relevant physical processes. Because they do not 
reveal the source of the nonlinearity, these studies have been of little value to the 
present work. 

C. Present Work 

In the present investigation we seek a simple theory which quantifies the 
effects of amplitude dependent resistivity. In particular, we seek to explain and 
rectify the discrepancies between Kuntz's theory and experiment. The theory is 
based on the following model of the porous material; the material is rigid, 
incompressible and homogeneous, and it has only two important properties, porosity 
Q and resistivity r. Knowledge of the microscopic structural details of the material 
is unnecessary. The porosity and resistivity are determined by direct measurement. 
The resistivity is shown to be the following function of velocity u: r=rjfr 2 usgn(u). 
The effects of iiydrodynamic nonlinearities, e.g., shock formation, are neglected, 
and sound propagation through the material is assumed to be isothermal. This model 
leads to a nonlinear wave equation, the solutions of which agree favorably with 
experimental observations. 

Two of Kuntz’s experiments liave been repeated (with minor modifi- 
cations): propagation of complex periodic signals, and dc flow resistivity measure- 
ments. Measurements were performed on batted Kevlar^^29 in the 400 to 6200 Hz 
frequency range for small-signal waves, and 500 to 1500 Hz for finite-amplitude 
waves up to 165 dB source sound pressure level (SPL). The porosity of the samples 
ranged from fl=0.94 to 0.9S. 


An outline of the current investigation is as follows. Chapter 2 begins 

with a section on material modeling. Equations of motion are then derived, which 

lead to a nonlinear wave equation and an infinite set of inhomogeneous Helmholtz 

equations. The solution of the Helmholtz equations and the predictions which can be 

made from them are given in Chapter 3. The chapter is divided into two sections: 

(1) propagation of low amplitude waves and (2) finite-amplitude effects. The s'.cond 

section i:. further divided into an analysis of the behavior of the fundamental and a 

study of the growth and decay of the harmonics. Chapter 4 is devoted to 

experiments. Discussions of the experimental apparatus and procedures are followed 

by a detailed comparison between theoretical predictions and measurements. 

Observed phenomena include excess attenuation of the fundamental, phase speed (of 

the fundamental) which decreases with amplitude, and a cubic distortion pattern in 

the harmonics. Hysteresis in the amplitude dependent phase speed seems to imply 

the presence of a nonlinear bistability. Chapter 5 is a summary and conclusion of 

the first four chapters of the study, which constitute the thesis work. Chapter 6 

was added after the thesis and presents an application of the theory to the acoustics 

of lined ducts. The vector wave equation is derived, and a frequency domain 

perturbation in the velocity potential is set up to solve the problems of reflection 

from a porous half-space and propagation in a lined duct. Appendix A gives a 

discussion of thermal effects on propagation in porous materials. A discussion of 

2 

the alternative use of a resistivity model r=rj+r^u is given in Appendix B. 
Appendix C contains the listing for and a short description of the computer program 
which calculates the numerical solutions. 


5 



CHAPTER 2 


THEORY AND MODELING 

A. Introduction 

The theories of sound propagation which pertain most directly to the 
present investigation were reviewed in Chapter 1. In Chapter 2 we develop a simple 
model for studying propagation of high intensity sound through air-filled bulk porous 
absorbers. Some of the assumptions made have been used previously by other 
authors, namely, 

1112 

1) the frame is assumed rigid, ’ 

2) dc flow resistivity can be substituted for the ac value through- 

*2,13 . 

out, ’ and 

3) isothermal sound speed is appropriate for all frequencies of 
interest.^ 

The fourth assumption, 

4) hydrodynamic nonlinearities are small relative to the resistive 
nonlinearity, 

is unique to the present work. 

In Section B the concepts of porosity and resistivity are defined. Assump- 
tions used in generating of the model are examined in detail. The continuity and 
momentum conservation equations, derived in Section C, are combined in Section C 
to produce a wave equation and the corresponding Helmholtz equation. 


Because the process has been assumed isothermal, the energy equation is not needed 
in the derivation of the wave equation and is therefore not discussed here. 
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B. Modeling 


Porous materials in general have geometrically disorganized structures; 
the frame elements ere arranged helter-skelter in near-random fashion. This fact 
makes detailed ana.ysis of the hydrodynamic field virtually impossible. We 
therefore prefer to take a macroscopic viewpoint in which the summed effect of 
local flows determines the material properties. 

1. Porosity 

The porosity is defined as the fraction of the material volume which is 
filled with air. We have devised an abstract function D(x,y,z) to describe the 
structure of a porous material; the function is merely a convenient mathematical 
construct and need not be directly measurable. The value of D is 0 at a point 
occupied by a structural element and 1 otherwise. It therefore follows that the 
F>orosity is simply the value of D averaged over a volume substaritially larger than 
the mean pore size; 

,y,z) dx dy dz. (2.1) 

The sample is assumed to be homogeneous and the porosity stationa'^y 
throughout the sample. In practice, the porosity is determined by empirical means; 
the known densities of a sample p^, the frame pj, and the air trapped within p^ 
combine in the relation 

^ =(Pf - P^)/(Pf - P^) . ( 2 . 2 ) 

While the porosity defines the fraction of volume occupied by air, the 
percent open area Q defines the unobstructed fraction of a representative plane 
through which the flow passes. First, the local percent area fl for a plane located 

3 

at z=Zq is 


8 


The average of the local value over all possible planes is the quantity <sj > , 

a 


Note that Eq. (2.4) reduces to Eq. (2.1), so that in general 

n • 

The percent open area and the porosity must be equal because the aggregate open 
area of all infinitesimally thin slices defines the total open volume. This result is 
obtained regardless of the material geometry. 

2. Resistivity 

Viscous drag causes the material to resist the passage of an acoustic 
wave. The resistance per unit length is called the resistivity r, and is defined for a 
one-dimensional flow as follows. 

dp/dx = -ru , (2.5) 

where dp/dx is the pressure gradient parallel to the flow and u the flow 
velocity. 

Since in general the details of the flow cannot be analyzed, the resistivity 

cannot be calculated explicitly and must therefore be measured directly. A dc flow 

test is the most straightforward way of meeisuring the resistivity. The fact that the 

dc flow resistivity also applies to oscillatory flows is not intuitively obvious, but 

2 

Zorumski and Parrot have shown that the ac and dc resistivities are indeed 
equivalent over a wide range of particle velocities. It is therefore possible to 
determine the ac resistivity of the medium by measuring the dc resistivity. 


Measured data indicate that the resistivity depends on the flow velocity; a 
mathematical expression of the resistivity is sought which is appropriate for both ac 
and dc flows. Static flow resistivity data^ (see Fig. 2-1) seem at first glance to 
indicate that the resistivity is of the form 


in which case 


r = rj + r^ u , 


( 2 . 6 ) 


-dp/dx = rjU + r 2 ^ 



(2.7) 


Closer examination, however, reveals the implausibility of using Eq. (2.7) for an ac 
signal: it suggests that the nonlinear contribution to the pressure gradient is 
positive regardless of the instantaneous direction of flow. If the discrepancy is to 
be remedied, the pressure gradient must not contain any even functions of u. The 
simplest improvement is to affix sgn(u) to the r 2 term in Eq. (2.7). The pressure 
gradient then becomes 

2 

-p^ itjU + r2U sgn(u) , (2.8) 

from which it follows that the resistivity has the form 


r = r j + r 2 u sgn(u) , (2.9) 

where u sgn(u) is the instantaneous speed of the flow. The coefficient 
r j is referred to as the linear resistivity because it plays a large part in determining 
the small-signal propagation effects. Similarly, the nonlinear resistivity is so called 
because it causes nonlinear behavior. Equation (2.9) is still compatible with the 


t 


A detailed discussion of the data is given in Chapter 4. 


The large number of trials 


at small flow velocities was taken to find ein accurate value of rj. 




dc flow resisUvity data. Although the addition of the sgn(u) term seems redundant 
for dc flows, its omission led Kuntz^ to inaccurate theoretical predictions of 
nonlinear behavior for acoustic signals. 

Other expressions exist which fit the data well and satisfy the pressure 
gradient requirement. One example is 

r = r 1 + r 3 u^ . (2.10) 

Equation (2.10), which we shall refer to as the "quadratic model" rather than the 
"sgn(u) model," is more attractive for analytical study than Eq. (2.9), but does not 
bring us any closer to closed form solutions of the wave equation. Furthermore, 
predicted behavior based on this model does not differ substantially from that 
predicted for the sgn(u) model (see again Fig. 2-1). Although we have carried out 
parallel theoretical developments using the quadratic model, the results are 
somewhat redundant and are not covered here in derail. The quadratic model can be 
used to obtain a perturbation solution. Analytical expressions for propagation of 
intense tones and numerical results for propagation of more complex finite- 
amplitude signals subject to quadratic resistivity are given in Appendix B. The 
sgn(u) model, however, receives the bulk of the attention in the remainder of this 
work. 

3. Fr ame Rigidity 

The material frame is set in motion by the sound field but under certain 
conditions is nearly stationary and can be considered rigid. Lambert and Zwikker 
and Kosten*® determined that forces coupling the frame and gas motion are 
insignificant if 

‘^l>'’l/Pes ’ <2.1 1) 
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where a>j is the decoupling frequency and is the effective structural density 
defined by Lambert. Interpretation of Eq. (2.11) is straightforward. The resistivity 
factor represents viscous coupling which drags the frame with the flow while the 
factor represents inertia which inhibits frame motion. For frequencies above U) ^ 
the inertia dominates and the frame tends to remain motionless. 

Kevlar*^29 was used extensively in the experimental phase of this inves- 
tigation. In the porosity range 0.94-0.98 tlie decoupling frequency is between 60 and 
85 Hz. The lowest fundamental frequency used in our experiments was considerably 
higher, namely )00 Hz. Thus, we have assumed the frame to be rigid. 

4. Constitutive Relation 

Heat transfer between the air and the frame causes the thermodynamic 
state of the entrapped air to depart from being adiabatic. We are primarily 
interested in the proper relationship of the pressure and density which, as the ideal 
gas law 

£ = £RT (2.12) 

:.:iows, is determined by the temperature fluctuations in the gas. The 
acoustic pressure p and acoustic density p are related to the total pressure £ and 
total density £ as £=Pq+P and £= p^ + p where p^ and p^ are ambient values. The 
total temperature is similarly expressed as T=T^+T and R is the gas constant. The 
frame, by vrtue of its considerable heat capacity, tends to remain at the ambient 
temperature T^ while the air temperature fluctuates with the acoustic signal. A 
thermal boundary layer forms around the frame. In the boundary layer the 
temperature fluctuations decay to near zero at the frame surfaces. The boundary 
layer thickness varies as the inverse square root of the frequency and, because of 


13 


the extensive surface area of the pores, may embrace significant fractions of ihe 
fluid. The mean temp>erature fluctuation in the pores determines the sound speed 

c = ydp/d p 

Although the phase velocity of a wave in a porous material is actually substantially 
less than the sound sp>eed, the following discussion is couched in terms of the sound 
speed because it is the most familiar name lor the quantity. Appendix A gives a 
detailed discussion of the effect on sound speed of heat conduction between the 
fibers and a laminar flow. For the present it is sufficient to summarize the results 
of Appendix A. 

The sound speed takes a frequency dependent value somewhere between 

the extremes of the adiabatic value c rt) To =vART and the iso*herrnal value 

o ' *^o o * o 

c.s^p^/ pQ= At high frequency, because the boundary layers around the fibers 
are relatively thin, most of the air is in an isentropic state. The isentropic value of 
the sound sp>eed is therefore appropriate. At low frequencies the thermal boundary 
layers fill the pores almost completely so that the air temperature is constrained to 
the ambient value T^. In this range the sound sp>eed value is roughly that for an 
isothermal gas. At sub-audio frequencies, however, the quasi-static temperature 
variation induced by the acoustic flow overcomes the heat capacity of the frame 
and drives the temperature of the material as a whole. Since the air and the frame 
are in thermal equilibrium, the adiabatic sound speed is once more appropriate. 
Figure 2-2 shows the approximate variation of sound speed with frequency for 
Kevlar^9 of porosity 0.94, 0.96, and 0.98. 

The sound speed calculation is based on knowledge of the temperature 
distribution around an individual fiber. This distribution in turn depends on the 
assumption that adjacent fibers do not interact thermally; that ’s, the boundary 


Q=0.98 

Q=0.96 

0=0.94 



0 
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layers must not overlap. If they do overlap the ar tual temperature distribution 
cannot be determined, but is expected to yield a sound speed that is (nore isothermal 
than Fig. 2-2 would indicate. The boundary layer interaction turns out to be signifi- 
cant at surprisingly high frequencies. The vertical liries which intersect the sound 
speed curves represent the frequencies below which the sound speed calculation 
cannot be trusted. A detailed discussion of how these frequencies are determined is 
given in Appendix A. 

The results shown in Fig. 2-2 are based on the assumption of laminar flow 
conditions. The high particle velocities of intense waves, however, render the flow 
highly turbulent. Carman*^ has determined that the anset of turbulent flow in 
porous materials occurs when the modified Reynold's number 




(2.13) 


reaches a value in excess of unity, where S=4(l-R)/d is the surface area of the fibers 


per unit volume, d the fiber diameter, and i> is the kinematic viscosity. For 
Kevlar 29 in porosities ranging from 0.98 to 0.94, the velocity at which turbulent 
flow begins is 0.1 to C m/sec, respectively. The flow velocities used in our 


experiments often exceed these values by as much as an order of magnitude: the 
assumption of laminar flow is clearly not applicable in these cases. Turbulent 


mixing of the fluid has a significant effect on the sound speed because it greatly 
improves the heat transfer between the air and the fibers. We therefore expect 
conditiens to be significantly closer to isothermal than would be predicted from a 
laminar flow model. Since we are primarily interested in me’isurements at high 
intensity in frequency ranges where the boundary layer interaction is significant, we 
are justified in assuming the isothermal sound speed. 
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5. Relative Importance of Hydrodynamic No nli nearity 


Hydrodynamic nonlinear effects are extremely important for finite- 

amplitude progressive waves in an open medium, but their influence on wave 

propagation in porous materials is small. Many fine discussions of the physics of 

18 1 9 20 2 1 

intense sound propagation in air already exist; » * * vve shall cover only those 

points pertinent to the present discussion. An intense wave in an open medium 
distorts because the propagation speed dx/dt varies over the waveform, 

dx/dt = + ^u(X ,t) , 

whe'^e ^ =(Y + 1)/2 and y is the ratio of specific heats. The waveform distortion, 
which we shall call hydrodynamic distortion, is cumulative. An initially sinusoidal 
wave tends to distort into a sawtooth wave, i.e., a periodic sequence of shock waves. 
The shock formation distance 

X = 1 //3ek^ , 

where e-’i /c is the dimensionless wave amplitude and k =27rf /c is the wave 
number, serves as a measure of the rate of distortion. x is large the distortion 
builds up very slowly; the reverse is true if x is small. The shock formation distance 
in our experinr.ents ranged from about 1.5 m to 5.5 km. These values should be 
compared with the characteristic length for the absorption, 

X = [/a. 

a 

where a is the small-signal attenuation coefficient. Since values of x^ for our 
experiments were between 2 and 9 cm, i.e., x»x , the wave was completely 

cl 

absorbed before hydrodynamic distortion had a chance tc develop. Furthermore, the 
effects of the nonlinear resistivity turn out to be much more severe than those for 
hydrodynamic distortion (see Chapter 3). Our experimental concitions are expected 
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to be typical of those encountered in practical uses of Kevlar as a high intensity 
sound absorber. We therefore assume the hydrodynamic nonlinearity to be negligible 
relative to the resistive nonlinearity and neglect any second and higher order terms 
not pertaining to the resistivity in the following derivation of the equations of 
motion. 

6. Structure Factor 

A true one-dimensional flow in a porous material exists only as a 
convenient fiction. In reality the flow is deflected by the frame in directions 
perpendicular to its net direction of travel. Some of the force applied in the x 
direction to a fluid particle is therefore wasted in the process of accelerating the 
fluid arcund obstacles. This inefficiency causes the fluid to appear more dense than 
it actually is. The effects of increased inertia are most notable at high frequency, 
where they cause the phase speed and attenuation to be reduced. 

Most studies of porous materials include factors to account for the extra 

inertia. They are called by many names; structure factor, apparent density, 

tcrtuosity, and others. Some of the theoretical discussions of the underlying physics 

10 15 22 

are quite elaborate. ’ ’ Direct measurement of the structure factor is, how- 

ever, not possible nor can one determine its value theoretically //ichout considering 
the microscopic structural details of the material. Furthermore, the dependence on 
amplitude of the structure factor is not known. It seemed better to this author to 
do without such an enigmatic parameter, lest in our desire to assign it a numerical 
value we reduce it to a "fudge factor". We therefore choose to overlook the 
structure factor in formulating the equations of motion, physical significance 
notwithstanding. 
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C. Equations of Motion 


The foregoing discussion of the physics of acoustic flow in a porous 
material demonstrates that the resistivity effects far outweigh those of frame 
motion, frequency dependent sound speed, and hydrodynamic nonlinearity. The 
latter effects are therefore neglected so that needless complication is eliminated. 
The problem is thus reduced to a rather simple one: plane wave propagation through 
a "fluid" with one easily definable prof>erty, namely amplitude-dependent resistivity. 
The purpose of this section is to convert this physical understanding into equations 
of motion. 

A number of general comments and definitions must be made at the 
outset. We are presently concerned with only one-dimensional aggregate flows. 
Local perpendicular flow components are assumed to sum to zero over a given plane. 
We define the control volume, depicted in Fig. 2-3, to be a differential section of 
length Ax and cross-sectional area A. The material within the volume is 
homogeneous and has the same porosity throughout. Flow enters the volume at the 
left side through the open area fiA and exits through the same area to the right. 
Pressures exerted on the sides of the control volume act on the entire cross- 
sectional area A, not just the open area. A resistive force due to vise js drag is 
presumed to act throughout the volume, not just in the pores. The velor density, 
and pressure are mean values across a given plane. The velocity and density are 
averaged only within the pores. 

1. Continuity Equation 

The equation of continuity mathematically expresses the balance between 
the rates of mass influx, efflux, and accumulation. At some instant in time the flow 
carries a certain amount of mass into the control volume through the open area at a 
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Figure 2-3 

Stationary control volume for determination of the 
conservation of mass relation for porous materials 
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rate (pu)(An)[^. Mass exits the volume at the rate ^x' expect 

the mass to be conserved, the difference between these two rates represents a rate 
of mass accumulation, p^HAAx. The continuity equation may therefore be written 
as 

p^flA AX = puflA|^ - puliAj^^^j^ . (2.14) 

The differential form of Eq. (2.14) is found by dividing by SlA Ax and taking the limit 
as Ax goes to zero; 

Pt+(pu)^ = 0 . (2.15) 

The porosity factor does not appear in this equation because the volume fraction in 
which the mass accumulates is equivalent to the area fraction through which the 
flow passes.^ 

Nonlinear terms not pertaining to the resistivity are neglected on the 
basis of the discussion in Section A. 5 . Furthermore, we assume that there is no dc 
flow through the sample. We can therefore reduce Eq. (2.15) to the familiar 
linearized continuity equation 

Pt + Pq\ = 0 • ^2.16) 

2. Momentum Equation 

The derivation of the momentum equation proceeds along similar lines. 
An identical control volume (depicted in Fig. 2-4) is established, but this time the 
momentum of the volume is considered. The rate of momentum accumulation 
(pu)^J2AAx is balanced by the net momentum influx to the volume and body and 
surface forces acting on the control volume. Momentum enters the volume at the 

^ Particle velocity at or within a rigid obstacle is, of course, zero. 
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Figure 2-4 

Stationary control volume for determinati n of the 
conservation of momentum relation for porous materials 
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left at the rate (pu)uBA|^ and exits to the right at the rate 
exerted on the sides of the volume result in surface forces £A|^ at the left side and 
gA|x^^x right. A body force due to viscous drag ruAAx acts throughout the 

volume in opposition to the flow. Although the form of r(u) is already known, it has 
been left in its general form for the moment. 

The mathematical statement of conservation of momentum is 

(pu)^fiAAx = (pu)u Af2|j^-(pu)u x-£A| AAx . (2.I7) 

This equation is divided by A Ax and take the limit as Ax goes to zero to find the 
differential form 

n(pu)^ +fl(puu)^ + + ru = 0 . (2.18) 

Again dropping hydrodynamic nonlinear terms, we obtain 

j]p u^ p + ru = 0 . (2.19) 

“^o t *^x 

The porosity factor appears in Eq. (2.19) because the flow enters the 
control volume through an area determined by the porosity while the pressure and 
body force act on the entire surface and volume, respectively. 

3. Pressure-Density Relation 

We recall from Section B.4 that the gas within the pores is constrained to 
the ambient temperature value T^ because of rapid convective heat transfer 
between the air and the frame. Equation (2.12) therefore reduces to 

p = pRT^ , (2.20) 

where c.^=jRT^ is the isothermal sound speed. 
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D. Wave Equaticxn 


In the previous sections the continuity, momentum, and constitutive 

equations appropriate for high intensity sound propagation through rigid bulk porous 

materials were derived. The wave equation is found by combining these three 

2 

relations. We first apply the differential operator c. 9/9 x to the continuity 
equation and 8/at to the momentum equation. Density and pressure terms cancel 


via Eq. (2.20), and the wave equations is found to be 


flu^^ - d?u + 

tt 1 XX 


(ru). 


( 2 . 21 ) 


We adopt the sgn(u) model for the nonlinear component of r. Substitution 
of Eq. (2.9) in Eq. (2.21) yields 


2 W „ 

”“tt - 'i“xx T “t * i 

o '^o 


( 2 . 22 ) 


It turns out that solutions of Eq. (2.22) are best sought in the frequency 
domain. If the nonlinear term is ignored for the moment, we see that the behavior 
implied by Eq. (2.22) is a combination of lossless plane wave propagation 


and a diffusion process 


Qu.. - c. u = 0, 

tt 1 XX ’ 


'"l 2 

— u^ - c/ u =0. 

t 1 XX 


For most cases of practical interest the diffusive behavior is dominant; that is, the 
resistance term TjU^/p^ is much larger than the flu^^ term. The diffusion equation 
has wavelike solutions, but they are characterized by severe frequency dependent 
attenuation and dispersion. In particular, no general time domain solutions like 


g(x>ct) or f(x-ct) for incoming or outgoing waves, respectively, can be found. 
Functions corresponding to incoming and outgoing waves are, however, identifiable 
in the frequency domain. It is therefore advantageous to transform Eq. (2.22) into 
its frequency domain counterpart, the Helmholtz equation. 

Before tranforming Eq. (2.22) we find it convenient to rewrite the 
equation in nondimensional form, vis., 

where 

v = u/Cj , ’ 

T = o;^t , = "2V^^o‘^o ’ 

X=w^x/c , 

and is a characteristic frequency, such as the fundamental frequency of the 
incident acoustic wave. 

2. Helmholtz Equation 

We now transform Eq. (2.23) to the frequency domain. The transfor- 

23 24 

mation is done in a manner popularized by Fenlon and Korpel. Time and 
distance dependences of the velocity variable v(x,t) are separated by expanding v as 
a doubly infinite series of exponentials in time, with range dependent amplitude 

•o 

v(X,t) = V^(X)exp(jnr) . (2.24) 

n=-«o 


The index n represents the nth harmonic of the signal. 
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Since v is a real quantity, we require that V ^ and be complex conjugates. The 
series approach is preferred to an integral transform because our incident signals 
are monochromatic or at least periodic. The nonlinear term contains the factor 

V sgn(v), here abbreviated as w(X,t), which represents the speed (not velocity) of the 
acoustic flow 


4* «o 

w(X,t)* vsgnv = ^ W^(X)exp(jnT) , (2.25) 

n = -«o 

and W ^=W^*. The speed w corresponds to a full-wave rectified version of the 
signal, and the spiectrum corresponds to the spiectrum of the velocity-dependent 
resistivity. 

Substitution of Eqs. (2.24) and (2.25) in Eq. (2.23) leads to the in- 
homogeneous Helmholtz equation. 


V + q V = jnSRO V W 
n n ' 2 / ^ p n-p 

P 


(2.26) 


where the double prime denotes two differentiations with respect to the argument 


X, and q^=n./fl(l-jR j/n) is the complex dimensionless wave number, which may be 
written 


q = q - I q . 
^n ^nr ’ ^ni 


(2.27) 


The real and imaginary parts are 
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(2.28) 

(2.29) 

(2.30) 

(2.31) 


3US Helmholtz equation is 
of Equation (2.26) are the 


CHAPTER 3 


SOLUTIONS 


A. Introduction 

Chapter 3 is devoted to solutions of the inhomogeneous Helmholtz 
equation for progressive waves. Solutions for low intensity linear behavior are 
investigated in Section B. Examination of these solutions gives a basic under- 
standing of the effects of resistivity on small-signal waves. An introduction to the 
theoretical problem of finite-amplitude waves in porous media is offered in 
Section C. An approximate first integral of Eq. (2.26) is developed that is the basis 
of the solutions in subsequent sections. In Section D, we investigate propagation of 
a tone subject to nonlinear resistivity effects. It is assumed here that the harmonic 
distortion products, thougli inevitable, are weak enough to be neglected. 
Concentratif'g, then, on the fundamental, we obtain approximate analytical solutions 
for pressure, particle velocity, phase velocity, and impedance magnitude and phase 
angle. Section E contains a discussion of the approximate solution for high intensity 
periodic waves in which the harmonics are taken into account. The multiplicity of 
harmonic interactions and interlocking solutions, as defined by the convolution term 
in Eq. (2.26), makes it necessary to seek a numerical solution. Harmonic 
propagation curves, i.e., amplitude versus distance for the various harmonic 
components, are generated. The numerical solution gives insight into the evolution 
of harmonic distortion products and in addition provides predictions against which 
experimental data may be compared. 
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A series of general expressions are now given in preparation for the 
analytical solutions discussed in Sections B and D. These definitions simplify the 
analysis and clarify the effects of resistivity in the linear and nonlinear regimes. 

1. Velocity and Pressure 

The particle velocity vfX,T) can be expressed as an infinite sum of 
sinusoids with range dependent phases and amplitudes A^(X); 

nr-0y^(X)^ , (3.1) 

where A (x)and o (x) are real and are related to V (defined in Eq. (2.24)) by 

n'"/ ’-yn n 



Use of the A^ notation allows the amplitude and phase effects to be easily separated 
so that they may be considered independently. 

The nondimensionalized pressure signal is defined as 


B(X,t)-p ^ ^ 

p(X,T) = = 2 ^ ^n^^^ exp(jnr) 

p/. n=.« 


A 

where P 


-n 



* and, alternatively, 


(3.3) 


CO 

p(X,T)= s n^(X)sin(nT-0p^(X)) 


(3.4) 


where n and <b 
“n ’^pn 


are real, and 
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(3.5) 



Hn(X) 

~W~ 



2. Phase Velocity 

A general expression for phase velocity, applicable in both the linear and 
nonlinear cases, may now be derived using the notation of Eqs. (3.1) and (3.4). The 
phase speed is defined as the rate of travel of a constant phase plane, that is, a 
plane for which the quantity nr- 0 y^(X) remains constant. Upon differen*iating this 
expression we find that 

ndr-o' dX=0 , 

'' vn ' 

where the prime again denotes differentiation with respect to the argument X . The 
phase velocity of the nth component of the particle vei>x:iiy wave is therefore 


C 


= ^ 
vn dr 


const. 

phase 


vn 


(3.6) 


The phase velocity for the pressure signal is similarly defined as the rate of travel 
of a plane for which nr- 0 pp(X) stays constant. 


r = 
^pn Hr 


const. 

phase 


’pn 


(3.7) 


Two phase velocities are defined because the impedance phase angle is amplitude 
dependent at high intensity and therefore the pressure and particle velocity signals 
travel at different speeds. 
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3. General Definition of Characteristic Impedance 

The characteristic impedance is defined in general as the ratio of the 
pressure to the prrticle velocity for a given harmonic component; 


7 - " 




(3.8) 


If is written in polar form, 


z„ = 




(3.9a) 


it follows that 


n„(X) 

l^nl ' 'VTxT 

n 


(3.9b) 


and 

j /-Im(Z A 

V®vn • 

A general relation for the imp>edance can be found by using the continuity 
equation which, although linear in form, provides a relation appropriate for all 
amplitudes because the hydrodynamic nonlinear terms have been neglected. We 
first replace in Eq. (2.16) by P^/Cj and t' en nondimensionalize the result to find 

Pr*“x=0 . 

Equations (3.1) and (3.3) are here substituted into this equation to give the 
dimensionless continuity equation in terms of the harmonic components. 
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(3.10) 


inP V =0 
‘ n n 


We first divide Eq. (3.10) by V^, 




and then substitute for V and S' using Eq. (3.2) to solve for Z, 

n n n 


Z 


n 


n 




(3.11) 


It is interesting to note that the rate of change of phase of the wave represents the 
real impedance while changes in the wave amplitude represent the reactive part. 

B. Low Intensity 

Attenuation and phase velocity are controlled at low intensity by the 
dimensionless linear resistivity coefficient Rj. The Helmholtz equation takes the 
linear form 

V"+q^V =0 (3.12) 

n ^n n 

for small signals. This second order differential equation may be factored into two 
first order relations by use of operational notation. Equation (3.12) is first written 
as 

or 

(D^.iqn)(D^.jqn)Vn = 0 , 

where represents differentiation with respect to x • The "roots" of this factored 
equation are 
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(3.13) 


aiid 


V jq V =0 
n n 


V’ - jq V = 0 . 

n n 


(3.14) 


Equation (3.13) is the relation for the forward traveling wave, Eq. (3.14) for a 
backward traveling wave. We now concentrate on Eq. (3.13), which may be 
interpreted as a first integral of the Helmholtz equation for forward traveling 
signals. 

Amplitude and phase information may be obtained by substituting Eq. (3.0 
into Eq. (3.13); one finds that 


A' ■*. j(q - 0 * )A = 0 . 

n ^ ^n '"vn n 


(3.15) 


The real and imaginary parts of this equation give relations for the amplitude and 
phase, respectively, vis.. 


and 


A’ + q .A =0 
n mi n 


q - o' =0 

^nr '^vn 


(3.16) 


(3.17) 


where q^^ and q^. are given by Eqs. (2.27) through (2.31). 

1. Particle Velocity Solution 

Solutions of differential equations (3.16) and (3.17) are quite simple. The 
rauige dependent amplitude is 


A^(X) = A^(0) exp(-q^.X) , (3.18) 

where A^(0) is the value of A^(X) at the effective origin of the signal. Thus q^j is 
the decay constant. The phase function is 
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(3.19) 


0 U) = q X ^ 0 (0) , 

'^vn ^nr'' ’ 

where 0yp,(O) ‘s the value of 0yp,(X) at the effective origin and we see that is the 
propagation wave number. The general solution for the particle velocity for small- 
signal progressive waves is therefore 

•O 

v(X,t) = A^(0) exp(-q^.X) sin^nT-q^^X- 0 y^(O)^ . (3.20) 

One effect of the resistivity is to attenuate the wave as it travels. The 
nondimensional rate of attenuation varies according to Eq. (2.29), with the dimen- 
sionless grouping Rj defined in Eq. (2.23). Figure 3-1 is a graph of Eq. (2.29) where 
for simplicity the value 0 = 1 heis been used. It is seen that the attenuation increases 
with the dimensional resistivity rj and frequency. The dependence of q^^ on 
porosity is complicated because r^ is itself a function of porosity. The range 
0<Rj<50 is sufficient to cover the audio spectrum for most materials. At high 
frequency (Rj«l) the attenuation, reaches an asymptotic value 

Q(=r j/2 p^Cj, which is independent of frequency. The low frequency asymptotic 
result is 


a 



I 


which is pioportional to </f just as in the case of boundary layer dominated 
attenuation. The high and low frequency results are not very useful, however, 
because the isothermal sound speed is inappropriate at very high or very low 
frequencies, nor is the frame rigid at low frequencies. 
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Another effect of the resistivity is to reduce the phase speed of the wave. 


The phase speed is evaluated by combining Eqs. (3.6) and (3.17) to obtain 


'vn ■ 0 ' 


(3.21) 


vn 


'nr 


Resistivity causes the phase velocity to be substantial!' less at low frequencies than 
the sound speed. The dependence of on Rj can be seen in Fig 3-2, where once 
again a value of n = l has been used for simplicity. The true phase speed value is 
c = The high and low frequency asymptotes, subject to frequency constraints 

on the theory, are the isothermal sound speed at high frequencies and 


c 



at low frequencies. The low frequency phase speed and attenuation have JT 
dependence which is indicative of a diffusion process. 

2. Small-Signal Characteristic Impedance 

The characteristic impedance is strongly affected by resistivity. 
Equations (3.16) and (3.17) are substituted into the general iinpedance relation, 
Eq. (3.11), to give 


Z 


n 




(3.22) 


At low frequencies sound propagation is reduced to a diffusion process, charac- 
terized by an impedance magnitude of s/7j/p^w and a phase angle of 45°. At high 
frequency the magnitude of the impedance approaches and reactance drops to 
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zero. Figure 3-3 illustrates the dependence of the impedence magnitude on Rj and 
Fig. 3-4 does so for the impedance phase angle. 

3. Small-Signal Pressure Wave Solution 

The solution for pressure waves is found by simply taking the product of 
the particle velocity and impedance. Since is independent of amplitude at low 
intensity, the pressure has a similar solution 

■O 

yx) exp(-q^.X) sin(nT-q^^X- 0 p^(O)) . (3.23) 


The phase velocity of the pressure wave and that of the particle velocity wave are 
therefore identical: 


P _ _n n_ 

pn " rt' ■ q 
'"pn ^nr 

The attenuation constant of the pressure signal is the same as well, 
expected result. 


(3.24) 
This is the 


4. Summary 

The effects of resistivity on low intensity acoustic signals are summarized 
as follows: the linear resistivity coefficient is the most important parameter 
controlling the flow of acoustic energy, since it 

1) causes attenuation, 

2) causes dispersion, and an overall decrease in the component phase 
speeds, and 

3) increases the impedance magnitude, and introduces a constant phase 
angle between the pressure and particle velocity signals. 
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Figure 3-3 

Dimensionless impedance magnitude versus 
dimensionless resistivity coefficient R, 





The method of separating the amplitude and phase information developed in this 
section is applied in Section D to the problem of a finite-amplitude tone. 

C. Introduction to the Nonlinear Problem 

The method of solution of Eq. (2.26) for finite-amplitude waves is dictated 
largely by the information available from the experiment. Without going into great 
detail, we can say that the experiment provides the harmonic pressure amplitudes 

A ) 

at several ranges for a progressive wavefield. While a single boundary condition 
per harmonic; is available, the inhomogeneous Helmholtz equation, Eq. (2.26), is 
second order and thus requires two boundary conditions per harmonic. The missing 
condition is the radiation condition, for which we have found no general analytical 
expression. Without this expression we can neither specify the second boundary 
condition exactly, nor can we integrate Eq. (2.26), in which case one boundary 
condition would be sufficient. Moreover, the missing boundary condition cannot be 
estimated accurately enough to prevent the numerical solution from becoming 
unstable. An approximate first integral is therefore developed which allows a stable 
solution for a progressive wavefield. Errors arise in the approximate solution but 
they are at least qualitatively predictable. 

For future reference we here introduce some special nomenclature. We 
are oftentimes concerned only with how the fundamental component of the wave 
propagates. The fundamental with higher harmonics neglected is, for simplicity, 
referreo to as a "finite-amplitude tone". When the higher harmonics are included in 
the analysis the signal is referred to as a "distorted sinusoid". 

1. Boundary Conditions 

The bcxindary conditions are taken to be the harmonic pressure amplitudes 
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of the signal at the microphone closest to the source, which we call the reference 
microphone. Although the transducer input is always sinus<Jidal, various nonlinear 
processes cause the high amplitude signal to distort before it reaches the first 
microphone position; the signal at the initial point therefore contains all harmonics 
of the original tone. The SPL of the hcirmonics are measured and from these the 
harmonic prer^ure amplitudes P^(0) are determined. These amplitudes serve as 

A 

starting points for the solutions PptX), to which we compare the data from 
microp>hones at other positions. 

In order to specify the solution of Eq. (2.26), it would be necessary to 
know the initial values (that is, the values at the initial point) of both V and V„'. 
However, only one of the boundary conditions can be synthesized from the pressure 
data; the value of V^' is linked to through the continuity equation, Eq. (3.10). The 
remaining boundary condition, however, is the value of at the reference 

A 

microphone. No exact general relation has yet been found between P^ and V^. The 
value of must therefore be estimated. We tried using various successive 
approximation methods to obtain the necessary estimates but found solutions of 
Eq. (2.26) to be so unstable that all attempted approximation schemes ended in 
failure. The reason for the extreme instability is that the Helmlwltz equation 
belongs to a class of differential equations called "stiff". In a stiff equation the real 
parts of the roots (eigenvalues) of the characteristic equation are widely separated, 
i.e., the values jq^.j are large. Attempted solutions based on one eigenvalue of the 
system are easily contaminated by errors (numerical and otherwise) which become 
forcing functions for the undesired solution. 

A better understanding of this problem may be gained by examining the 
homogeneous solution of Eq. (2.26), 

V^(X) = B^(0) exp(jq^X) + ®xp(-jq^X) . (3.25) 
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For a progressive forward-traveling wavefield B^(0) is expected to be identically 
zero. Unfortunately, we are unable to accurately specify the boundary conditions 
for finite amplitude waves, so that retains some residual value. Since the 

backward-traveling wave decays for decreasing X , it grows just as rapidly in the 
forward direction as the forward-traveling wave decays. Eventually a point is 
reached where the spurious backward-traveling wave solution obscures the desired 
solution. The larger the values of q^. the more rapidly the accuracy of the solution 
deteriorates. 

In order to achieve a stable solution, the second order system of equations 
must be reduced to first order by applying the condition that the wave field is 
progressive. An exact first integral of Eq. (2.26), similar in form to Eq. (3.13), has 
not been found because of the complexity of the convolution term. We therefore 
seek an approximate first integral. 

2. Approximate First Integral 

The approximate first integral isolates forward wave behavior within 
Eq. (2.26), but does so at the expense of sorr* ’ccuracy. Since we seek solutions lor 
the forward-traveling waves, we propose a solution of the form 

V^(X) = F^(X) exp(.jq^X) , (3.26) , 

where the deviation of F ( X ) from a constant value is due to nonlinear effects, and 
the exponential term represents small-signal effects for a forward-traveling wave. 
One may think of Fpj(^) as an apparent harmonic source amplitude which is allowed 
to vary with distance in response to nonlinear effects. Substitution into Eq. (2.26) 
gives 

''n" ■ = '"““z S''p*n-p 

P 
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Note that this substitution has not reduced the number of eigenvalues in the 
problem; to accomplish this feat an integration is necessary. The danger of 
contamination by spurious backward-traveling waves still remains. We have, 
however, shifted the eigenvalues of the linear problem to 0 and 2jq^, so that we can 
take advantage of the differing relative importance of the terms F^" and F^' for the 
two wave types. The following analysis is carried out as if there were no backward 
waves, and it is this selective vision which accomplishes the integration. For 
forward-traveling waves where the nonlinearity is moderate, the trcinsition between 
nonlinear and linear behavior is thought to be sufficiently gradual to warrant the 
assumption that F^" can be neglected relative to Removal of the F^” term 

results in a truncated version of Eq. (3.27), 


QR- 

F ’(X) = -n -: 5 -^ 'V' V W expMJ 
n 2q^ ^ P n-p ^ -'^n 


X) 


(3.28) 


In effect an approximate first integral of Eq. (2.26) has been obtained. 


V ' 
n 


jq V 
'^n n 



(3.29) 


The fact that Eq. (3.29) is the appropriate first integral is proved by noting that the 
relation l^n'I^^P^n^n'l satisfied by a backward- traveling wave. 

Equation (3.29) is therefore an approximate differential equation governing pro- 
gressive waves. 


3. Solution Errors 


In the process of rejecting the spurious solution, a certain amount of 
information regarding the spatial dependence of has been lost. Implicit in 




Eq. (3.29) is the supposition that the transition between nonlinear and linear 
behavior is grrdual. For very large amplitude waves, however, the transition is 
relat.vely rapid and F^" is not negligible. The errors of the approximate solution 
that result are called "overshooting" after their graphical manifestation. Non- 
linearly induced gains and losses for the fundamental and various harmonics (see 
Section E) of the signal are predicted too large. Since the exact solution is not 
known, there is no way to quantify the errors of the approximate solution. The 
errors will, however, become apparent when the solution is compared with measured 
data in Chapter 4. 

D. Propagation of a Single Tone 

We now consider propagation of a single tone of finite amplitude. The 
effects of the inevitable harmonic distortion components on the tone are neglected 
because the highe- harmonics are very weak relative to the fundamental. With the 
aid of the approximate first integral, we are able to find analytical expressions 
governing the amplitude and phase of the wave and the impedance of the medium. 

Analytical solutions of Eq. (3.29) are obtained for a tone 

v(X,t) = Aj(X) sin(T-0yj(X)| , (3.30) 

which has spectral components Vp(Aj/2j)exp(-j0^j(X)), and V j=Vj*. Since these 
are the only harmonic components assumed to exist, the convolution term of 
Eq. (3.29) is shortened considerably; 
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The spectrum of the full-wave rectivied version of an arbitrary signal 
is usually too complicated to analyze. In the pure tone case, however, it is 
straightforward. The rectified signal is defined as 


Aj(X) |sin( T- 0 ^j(X))| 
Application of the Fourier transform gives 




A.<X)exp(-jnT', .(X)) 

^ /sin(T')exp(-jnT*)dT' - j sin(T')exp(-jnT')dT' 


where t = t- 0 yj(X). Finally, evaluation of the integrals yields 


2A (X) 

W (X) = exp(-jn 0 j(X)) 

" TTd-n^) ^ ’ 


(3.32) 


for n=even and 0 otherwise. 

A differential equation in Aj and 0 ^j arises when we substitute Eq. (3.30) 
and (3.32) into Eq. (3.31), 


+ j(qi-0'vi)Ai 


- J2R 


2-fL 


q, 37 t • 


(3.33) 


As in the linear analysis (Section B), the real and imaginary parts of Eq. (3.33) give 
separate relations for amplitude and phase information. 


A' q,jA, .-Kq,^A, . 


(3.3«a) 


and 


^Ir ■ <®vl = ■’^^li'^1 ’ 


(3.34b) 
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where 


•rrv ^ 

K = (3.35) 

3ir^l + Rj 

is the coefficient of resistive nonlinearity. For the Kevlar samples used in our 
experiments the value of K was often 20 or greater. Since the hydrodynamic 
distortion effects would have been scaled by the dimensionless coefficient /?=1.2, we 
see that it is indeed correct to neglect them. 

The amplitude of the particle velocity wave is foun'-i by integrating 
Eq. (3.34a) directly, 


Aj(X) = 


Aj(0) exp(-qj.X) 

I + rJl-exp(-qjjX)j 


(3.36) 


where r=KAj(0)/Qj, which is directly analogous to the Gol'dberg number.^ 
Figure 3-5 gives plots of log(Aj) versus distance for various initial particle velocity 
amplitudes. Low intensity waves experience simple exponential decay. Hence, the 
corresponding curves are straight lines. High amplitude waves, however, undergo 
excess attenuation as a result of nonlinear resistivity. As a result, the propagation 
curves bend. For very high source amplitudes the pressure amplitude approaches an 
asymptote (the dotted line in Fig. 3-5); the amount of energy that can be 
transmitted to a given point cannot be exceeded regardless of the starting 


^ The Gol'dberg number, which was originally developed for distortion due to 
nonlinear effects in nonporous fluids, measures the relative importance of the 
nonlinear distortion effects to small-signal attenuation. 
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amplitude. The "saturation" limit is determined by the properties oi the medium 
and is independent of the initial amplitude. 

Although some have heretofore thought the excess attenuation in porous 

materials is due to conventional nonlinear effects in the fluid, we now see that the 
true cause is nonlinear resistivity. Equation (3.36) is of the same form as the 
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equation derived by Webster and Blackstock to describe excess attenuation and 
saturation of a nonlinearly propagating pure tone in free air. In their model, excess 
attenuation and saturation are a direct consequence of shock formation. Intense 
sound waves in porous media, however, do not form shock waves because the 
harmonics are rapidly attenuated and dispersed. The form of Eq. (3.36) was also 
obtained by Kuntz^ as a result of his empirical model of attenuation of a tone in a 
porous medium. The fact that Eq. (3.36) can be derived analytically proves that the 
basis of Kuntz's model was correct. 

The phase speed of the particle velocity wave is affected by the 
nonlinearity in a rather unexpected way. Substituting Eq. (3.10) into Eq. (3.6), we 
obtain 

The phase speed is seen to depend on the inverse of the propagation wave number qj^ 
at low intensity, but the amplitude dependence causes to decrease at high 
amplitude (see Fig. 3-5). A high intensity wave therefore starts by propagating 
slowly, but as the wave is attenuated the phase velocity increases, eventually 
reaching the linear value The resemblance of the characteristic shape of the 

curve in Fig. 3-6 and that in Fig. 3-2 should be noted. 

In Chapter U the relative phase between two points is measured as an 
indication of the phase speed, to which it is closely related. Equation (3.34b) can be 
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AS-85-568 


Figure 3-6 

Dimensionless normalized phase speed versus dimensionless 
particle velocity amplitude for fundamental component 




integrated to give the phase function 


<>vl(X) = q 


Ir 


X ♦ 




n (l + r (l-exp(-qj.x))) 




(3.38) 


The nonlinearly induced decrease in the phase speed is manifested here in the 
logarithmic term, which represents the increased relative phase at high intensity 
between two points. The full particle velocity solution is therefore 


v(X,t) = 


Aj(0) exp(-qj.X) 
1+ ^l-exp(-qj.x5j 


sin^T-qirX-Qiln (l+ r(l-exp(-qiiX)) - 


(3.39) 


The impedance of the wave is strongly amplitude dependent. The general 
impedance relation Eq. (3.10), combined with Eqs. (3.34a) and (3.34b), gives 

Zj = (qj^+KqjjAj) - j(qjj+Kq^^Aj) , (3.40a) 

or 


Zj = qj - jKAjq* 


(3.40b) 


The magnitude of the impedance can, after some manipulation, be shown to be 

1/2 


lz,| =^|q,|2(u(KA,)^).2R,KA,y 


(3.41a) 


which can be approximated by 

|Z||=|q,|(l.KA, ) 


(3.41b) 


when Rj is sufficiently greater than unity. As one may readily observe from 
Eq. (3.41b) and Fig 3-7, the impedance magnitude increases linearly with particle 
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Dimensionless normalized impedance magnitude versus dimensionless 
particle velocity amplitude for fundamental component 



velocity. The phase angle of the impedance is 
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zl 



/ qH>Kq, rAi\ 

Vq,r*Kq,iA|/ 
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Once again the phase information gives a surprising result. The impedance phase 
angle varies between a high intensity asymptote tan'^(l/Qj) and a low intensity 
limit tan'*(Qj). The material is relatively more reactive at high amplitude; as the 
wave travels and is attenuated, the reactance approaches its low intensity value. 
The amplitude-dependent impedance phase angle has interesting consequences for 
tlie pressure wave phase velocity, as we shall see shortly. A plot of Eq. (3.42) 
appears in Fig 3-8. 

The solution for finite-amplitude pressure waves is again simply the 
product of the particle velocity and inipedance, 

P(X,t) = AjV ) Zj(X) sin^T-0yj(X)- 0^j(X)j . (3.43) 

The pressure signal is subject to excess attenuation cind saturation just as the 
particle velocity is. But since the impedance magnitude is amplitude dependent, the 
functional dependence is slightly different. Although not apparent from Fig. 3-7, 
the impedance reaches a high amplitude limiting value which is determined by the 
particle velocity amplitude at saturation. The product of the saturation limits of 
the impedance and particle velocity as a function of x gives the high intensity 
pressure amulitude asymptote. Propagation curves for various initial pressure 
amplitudes are given in Fig. 3-9. 

The most surprising result of this section comes when we consider the 
phase velocity of the pressure wave which, from Eq. (3.7) and the definition of 
impedance phase angle, is 









(3.44) 


This phase speed is different from that for the particle velocity because at large 
amplitudes does not equal zero. The value of 0^j has previously been given, 
but the value of must be obtained by differentiating tlq. (3.42) with respect 
to X . It can be shown that the rate of change in the impedance ohase angle with 
distance is 


where 




-q|.Q|B exp(-qj.x) 


/ 7 

(3.45) 

fl - B exp(-qjjX)j + Qj 



r(i-Q.^) 

B = — 

1 + r 


The 0 ^'j term opposes the term. The pressure signal therefore travels slightly 
faster than the particle velocity signal. As the wave is attenuated, the pha*- speeds 
cf the two waves eventually converge to their mutual small-signal value In 

Fig. 3-10 the pressure phase velocity Cp^ is shown as a function of distance for a 
given source level. The corresponding curve for is provided for comparison. 
This phenomenon has been observed during the course of our experiments, and is 
reported in Chapter 4. 

Several general conclusions can be made regarding the effects of high 
intensity on propagation of a finite-amplitude single tone: 

1) The nonlinearity causes excess attenuation and, in the limit, 
saturation. 

2) The phase velocity is less than its small-signal value. 
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3) The impedance magnitude and reactance increase with amplitude. 

't The phase velocities of the pressure and particle velocity signals 
differ. The difference is due to the dependence of the reactance of 
the material on amplitude. 

Since the wave decays as it travels, the attenuation, impedance, and phase velocity 
eventually revert to their small-signal values. 

It is clear that small-signal resistivity effects are augmented at high 
amplitude, that is, the higher amplitude wave encounters a larger resistivity. The 
increased resistivity can be demonstrated by direct calculation. We consider a 
nondimensionalized version of Eq. (2.9), 

R = Rj + R 2 w(X,t) . (3.46) 

The average value of the resistivity over a cycle of the tone is simply 

<R> = Rj + R 2 W^(X) . (3.47) 

This result demonstrates a local increase in the mean resistivity because of the 
passage of the wave. 

The reader should note, however, that the resistivity actually fluctuates 
in time because of the acoustic signal, as defined in Eq. (3.46). Time-varying 
properties of the medium signify that new harmonic components are produced by the 
interaction of the signal and the fluctuating resistivity. In the next section we 
expand the analysis to include the harmonic distortion products as well. 

E. Numerical Solution for Initially Pure Tones and Distorted Sinusoids 

The purpose of the numerical solution is to generate theoretical pre- 
dictions for comparison with actual signals. We must therefore extend the solutions 
of the inhomogeneous Helmholtz equation to include the harmonics of the signal. In 
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particular, we are interested in the solutions for a distorted sinusoid like those used 
in our experiment. Closed form analytical solutions are not possible in the general 
case because the spectrum of the full wave rectified signal (W^) cannot be expressed 
as a direct function of the signal harmonics V^. In addition, the full solution 
consists of a large number of interdependent individual solutions, one for each 
harmonic in the signal. For these reasons the general case is best solved by 
computation. Nevertheless, qualitative predictions can be made by inspection of 
Eq. (3.29). 

A single initial value for each harmonic is required to solve Eq. (3.29), 
namely V^(0). The values of V^(0) are easily determined from the available data for 
P^(0) by means of ein approximate impedance relation that we now derive. The 
continuity equation, Eq. (3.10), relates and V^', ard the first integral Eq. (3.29) 
relates V^' and V^. These relations are combined to give an approximate impedance 
relation 


A 



(3.48) 


Since the distortion spectrum V W depenu’b on the signal spectrum V , we 

p P -P n 

determine V^(0) by successive approximation. Minute approximation errors are no 
longer a cause for concern, since the possibility of spurious wave solutions has been 
eliminated; the single initial value per harmonic is sufficient to fully specify a 
stable solution. 


The velocity dependent resistivity causes energy to be exchanged between 
the various harmonics of the signal. Two basic but loosely defined classes of 
harmonics arise: primaries and distortion products. Primaries are net suppliers of 
energy to other harmonics. They are characterized by the relation 
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(3.49) 


;fW I 

dx n' 



f 


which implies :i.at the attenuation rate exceeds that for a small-signal wave q^j. 
Other harmonics are net receivers of energy. They are called distortion products, 
and are characterized by the relation 





(3.50) 


Occasionally a distortion product receives such a rapid influx of energy that the 
harmonic amplitude initially increases with distance, 


S l''nl 


> 0 


(3.51) 


The harmonic amplitude does not increase indefinitely, however. As the primary 
decays the loss mechanisms soon overtake the energy influx. In this case the 
corresponding propagation curve takes on a characteristic rainbow shape. 

1. Propagation and Distortion of an Initially Pure Tone 

The case at hand is that of an initially pure tone, for which the rectified 
spectrum has been shown (see Eq. (3.32)) to consist entirely of even harmonics. In 
the case of a pure tone, the resistivity vat ies in time at even multiple frequencies of 
the fundament;'. As a result, odd harmonic distortion products are introduced into 
the signal. Despite their introduction, it turns out that the spectrum retains its 
even harmonic character. An initially pure tone is therefore expected to acquire 
only odd harmonic distortion components. This pattern is commonly referred tc as 
cubic distortion. 
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For the pure tone case, it is rather easy to recognize the primary and the 
distortion products. The differential equations gc ming the first four harmonics 
are 


V. 


jqiVi 


fiR- 


.) 


Vj(0) = 


10 


(3.52a) 


V '2 + jq2V2 = 0 ; V2(0) = 0 


(3.52b) 


V, 


j^3''3 = - 


■3nR 


-(VjW2+V*W^+V3W^+V5Wg+ ) ; V^(0) = 0 


(3.52c) 


and 


V’^ + jq^V^ = 0 ; V^(0) = 0 


(3.52d) 


The fundamental can be identified as a primary by tne term R 2 VjWq on the right- 
hand side of Eq. (3.52a). The equations for the even harmonics are homogeneous. 
Because the initial values for these equations are zero, the even harmonics stay at 
zero amplitude for all distances. The odd harmonics, however, are vigorously 
energized by the ''indamental as shown by the presence of forcing terms containing 
and Vj*. Because their initial amplitude is zero, the odd harmonics are expected 
to have the characteristic rainbow shape. Figure 3-11 is a graph of the propagation 
curves for an initially pure tone. It is not possible to compare these predictions with 
data because we were not able to produce a true pure tone at high intensity in a 
porous material. This example is, however, very useful for helping us understand the 
distortion of more complicated signals. 
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Predicted propagation curves for an 
initially pure finite-amplitude tone 





2. PropaRation of Distorted Sinusoids 


We now consider the case of propagation of distorted sinusoids like those 
used in our experiment. These signals are characterized by a strong fundamental 
accompanied by both odd and even harmonics which decrease in amplitude with 
increasing harmonic number. Although it is difficult to demonstrate analytically, it 
turns out that the addition of even harmonics to the signal has the effect of causing 
odd harmonics to appear in the rectified signal spectrum. We can make predictions 
of the harmonic propagation curve shapes by examining the differential equations 
governing the first four harmonics, 

V'l . jqiV, 

- -5^(V^W^+V|W2+V2W|+V5W3+V3W5+V5W^q- . . .) ; V^(0) = Vjq , (3.53a) 


VL + i 




2nR- 


■^(V2W^W*W^+ViWi+V|W3+V^W*+V*Wg+ . . .) ; V2(0) = V20 , (3.53b) 


-3RR. 


^3 * ^*^ 3^3 “ 
and 

V’. + 


<VjW 2+V|W^+V3W^+V5W^+ . . .) ; W^{0) = V3Q , (3.53c) 


-4HR, 




Again, only the most prominent forcing functions have been listed. 

The fundamental is the strongest component in the signal. It is subject to 
excess losses because of the augmented resistivity, but some of its lost energy is 
upshifted in frequency and strengthens the family of odd harmonics (the rest of the 
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energy is dissipated). The second harmonic is less strong than the fundamental but, 
because a negligible amount of energy reaches it from the fundamental, it too 
functions as a primary. Most of its energy is lost to augmented resistivity and in so 
doing it energizes the higher even harmonics. In fact, a second harmonic 
independently undergoing cubic distortion would energize the sixth, tenth, and 
fourteenth harmonics, and so on. It is therefore of interest to note that the fourth 
harmonic and its multiples are not members of the second harmonic cubic distortion 
product group. They are products of coupling between the primaries as represented 
by the forcing function '^ 2^2 relies heavily on the fundamental strength 

through the factor W 2 . Figure 3-12 shows the first five numerically calculated 
propagation curves for a distorted sinusoid. The initial harmonic levels are those of 
an actual measured signal. 


3. Cubic Nonlinearity 


All of the aforementioned distortion spectra could be explained equally 
well in terms of a cubic nonlinearity, i.e., the nonlineai’ity which results from use of 
the quadratic resistivity form 

2 

r = r j + r j u 

In this case, the corresponding form of Eq. (3.29) is 


V' + jq^V„ = 
n '^n n 


-nSlR 


V, 


(3.5*) 


Analysis of the convolution term reveals that a pure tone input produces odd 
harmonic distortion products. Indeed, the "cubic" distortion pattern was named 
after the effects of just such a nonlinearity. This pattern leads us to suspect that 
other results for cubic distortion are similar to those for v sgn(v) distortion. 
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The computer program for the numjrical solution is easily reconfigured to 

accommodate the cubic nonlinearity. The modified program was run starting with 

the same initial level*-, used to produce the curves in Fig. 3-12. Tne predicted 

behavior for cubic nonlinearity, shown in Fig. 3-13, is qualitatively very much like 

that for the v sgn(v) nonlinearity. This should come as no surprise, since the 

resistivity curves r(u) for the two models are designed to mimic the same set of 

resistivity data. Although the u sgn(u) model gave a better fit to the resistivity 

data, the quantitative results are not that much different here. Moreover, when 

compared to the measured acoustic data the differences between the two solutions 

are smaller them the approximation errors inherent in the solution (see Chapter 4 

and Appendix B). We are therefore led to believe that a cubic nonlinearity may be 

2 

substituted for the v sgn(v) nonlinearity. This approach will be of considerable 
usefulness for a planned perturbation solution of a high intensity sound field in a 
lined duct. 

Analytical solutions of Eq. (3.54) for the finite-amplitude tone case are 
very similar to their counterparts based on Eq. (3.29), and are given in Appendix B. 

We have acquired basic understanding of the harmonic distortion caused 
by nonlinear resistivity. The numerical solution, which is based on an approximate 
first integral of Eq. (2.26), provides quantitative predictions foi the cases of initially 
pure tones and distorted sinusoids. The fundamentals of the distorted sinusoids are 
shown to behave like the tones discussed in Section D, cmd the growth of harmonics 
is in a cubic distortion pattern. In Chapter 4 we compare the numerical solution 
results with measured data. Details of the computer program used to implement the 
solution are discussed in Appendix C. 




CHAPTER 4 


EXPERIMENTAL WORK 


A. Introduction 

Three experiments have been performed to measure the effects of a 
porous material on static flows and acoustic signals. Comparison of measured data 
to theoretical predictions provides an indication of the validity of the theory. We 
discuss the experiments in the following order;^ 

1) static measurement of flow resistance, 

2) propagation measurements at high and low intensity, 

3) phase speed measurements. 

In the first we measure the pressure drop across a length of bulk porous material for 
various dc flow rates, from which we determine the resistivity of the sample as a 
function of flow speed. The second experiment comprises measurements of the 
harmonic sound pressure levels at various ranges. In the third test we determine the 
phase speed from the relative phase of the fundamental between two points. The 
latter two experiments are referred to here as the acoustic tests. 

The dc data and the acoustic data serve different purposes. By fitting 
Eq. (2.9) to the dc flow resistivity data for a given sample, we determine the 
constants rj and r 2 ^ These constants are treated as properties of the medium and 
serve as input parameters for the numerical and analytical solutions. The 

experiments are intended to confirm or refute two of our key assumptions, namely: 


The actual chronological order is discussed in Section C of this chapter. 
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1) the dc and ac flow resistivities are the sa^ne, and 

2) the dominant nonlinear mechanism is a speed dependent resistivity 
of the form r=r j+r 2 usgn(u). 

The ac resistivity inferred from the small-signal attenuation data is compared with 
the dc resistivity to test the former assumption, and the qualitative agreement of 
the harmonic data and predictions is taken as an indication of the validity of the 
latter assumption. Data from the phase velocity test provides, for comparison with 
the prediction of Eq. (3.35), that a high intensity tone propagates more slowly than a 
small-signal wave of the same frequency. 

The static flow and acoustic tests previously performed by Kuntz used a 
separate apparatus for each test type. Hence the static flow and acoustic 
measuremer. s were performed on different samples. Since it is very difficult to 
create samples with identical properties, one is led to wonder about the com- 
patability of the two data sets. We have therefore modified Kuntz's traveling wave 
tube to accommodate both the acoustic and static flow experiments. Experimental 
results indicate that the key assumptions stated above are at least qualitatively 
valid. 

The remainder of Chapter 4 is divided into four parts. Section B contains 
a discussion of the equipment used in the experiments. Experimental procedure is 
detailed in Section C. Finally, the experimental data and theoretical predictions are 
compared in Section D. Section E is a summary of the experimental work and the 
conclusions drawn therefrom. 


70 


B. Experimental Apparatus 


This section is divided into discussions of the traveling wave tube, the dc 
flow teot apparatus, and the acoustic test apparatus. Since the traveling wave tube 
is common to all experiments, it is discussed first, followed by a discussion of the dc 
flow apparatus and of the acoustic test equipment. 

The traveling wave tube is the central piece of equipment for all the 
experiments. The sample fills the tube and is held in place by friction with the 
walls. The tube provides for measurement of acoustic and static flows in the 
sample. The length of the tube is 26 in., the inside cross-section is 3/4 in. on a side, 
and the wall thickness is 1/8 in. A screen of 0.008 in. steel piano wire restrains the 
sample just inside v'hat we call the downstream end of the tube. A series of H in. 
diam. holes are drilled in the top side of the tube to serve as microphone ports. 
Engineering drawings of the tube are presented in Fig. 4-1 for the convenience of 
the reader (reproduced from Kuntz's dissertation by his kind permission). 

2. dc Flow Apparatus 

The dc flow apparatus measures the pressure drop across a porous sample 
as a function of flow speed. The traveling wave tube was originally designed for 
acoustical measurements, not for high pressure dc flow tests, so the tube had to be 
modified to accept a static flow input and to prevent leakage from the microphone 
holes. An adapter, fabricated from a block of hexagonal aluminum stock, is affixed 
at the upstream end of the tube to provide attachment points for the incoming flow 
and for the pressure gages. The adapter is clamped to the tube and a gasket 
guarantees a tight seal. A strip of rubber gasket material backed with aluminum is 
clamped on top of the tube to seal the microphone holes. These modifications were 
quite effective; at no time was the tube observed to leak. 






The apparatus is best understood if the flow is followed from source to 
outlet (see thie block diagram in Fig. 4-2). A regulated 100 psi air supply drives the 
flow. The flow rate in standard cubic feet p>er hour (SCFH) is measured by one of 
a collection of flowmeters which cover various ranges. The flow then enters the 
tube at the adapter, passes through the sample, and exhausts into the atmosphere at 
the left end of the tube. The gage pressure in inches of water (in. H 2 O) at the 
adapter block is measured by a series of pressure gages. Since the gage and tube 
snare the same outlet pressure, the gage reading indicates the pressure drop across 
the sample Ap. The resistivity is 


r 

Fr'- 


where L is the length of the sample and A is the cross-sectional area of the tube. 


3. ac Test Apparatus 

The ac test apparatus permits measurement of harmonic SPLs and 
relative phase of the fundamental at six down range positions. Six microphones 
were used in the experiment so that they would not have to be moved repeatedly. 
The only drawback to this afjproach is that the signals from the microphones must be 
adjusted so that the microphones respond equally to a given input. The microphone 
nearest the source is called tiie reference microphone; the harmonic spectrum 
measured at this position provides the boundary conditions for the numerical 
calculation. Phase angle measurements are referred to the signal ai this micro- 
phone as well. The ac apparatus includes some rather elaborate electronics; the 
reader is encouraged to make use of the block diagram provided in Fig. 4-3. Four 
distinct sections of the apparatus are identified; 
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FIGURE 4-3 
BLOCK DIAGRAM OF 
ACOUSTIC MEASUREMENT APPARATUS 
























1 ) 

2 ) 


the source section, 
the acoustical section, 

3) the measurement section, and 

4) the control elements. 

They are discussed in this order. 

The source section begins with the oscillator and ends with the com- 
pression driver.^ The oscillator produces a continuous pure tone which is appro- 
priately amplified or attenuated by the voltage controlled amplifier (VGA) section 
of the compressor. Thie compressor is a control element and is discussed below. 
Since signal distortion is intioduced by the VGA, a high-Q bandpass filter tuned to 
the fundamental signal frequency is inserted in the circuit to remove unwanted 
harmonics. It is desirable to have as pure a spectrum as possible at the reference 
microphone because the harmonic distortion effects are more apparent. A decade 
attenuator compensates for the gain of the filter. The source signal then passes to 
the DuKane 200W amplifier and on to the acoustic source, a 3BL 375H compression 
driver with titanium diaphragm. The power applied to the driver was kept below 
25 W continuous rms which, for 8 52 nominal impedance, corresponds to a current 
limit of 1.768 A as measured by the Keithley VOM. The oscilloscope was used to 
visually warn the experimenter in the event of signal fluctuations or clipping. 

The acoustic section of the experiment includes the tube, all associated 
hardware, and the sample itself. The empty tube has a first cross- mode cut-on at 
about 9 kHz; frequencies used here are low enough that inadvertent m<;dai 
excitation is not expected to be a factor in the measurements. An adapter connects 


^A compression driver is usually used to drive an acoustic horn. 


the driver to the tube by means of a tapered bore which smooths the cross-sectional 
area difference. The microphones are placed in holders at the reference (0 cm), 2 , 
5, 8, 11 and 20 cm positions. In addition, the microphone holders grip the tube 
tightly. Thus, tube wall mass loading is provided, which reduces vibration of the 
holders. A narrow air gaj ' left between the microphone face and the sample to 
avoid noise caused by direct contact of the fibers with the microphone face. 
Finally, the porous sample serves as its own acoustic termination, since it provides 
small-signal attenuations in the hundreds of dB/m. The sample extends well beyond 
the 20 cm microphone position; no reflected waves are expected to be present 
within the sample near the microphones. 

Although the source signal is originally sinusoidal, the signal received at 
the reference microphone is somewhat distorted at high intensity. Harmonics are 
added to the signal by a variety of nonlinear mechanisms, including transducer 
nonlinearity, propagation through the air, transmission through the material/air 
interface, and propagation through the material. 

The measurement signal paths begin at each of the six microphones. Five 
Bruel and Kjaer model 4136 and one model 4135 microphone were used. The signals 
from the individual microphones pass to a rack of instrumentation amplifiers where 
the dc power supply voltages are removed and coarse microphone calibration is 
done. Fine calibration is performed with a separate set of potentiometers so that 
the microphones respond uniformly to the calibrated signal from a B&K pistonphone. 
A six-position rotary switch is used to select the output of one microphone channel 
at a time. The harmonic spectrum of the signal from the selected microphone is 
measured with either a B&K 2010 heterodvne analyzer or an HP 3580A spectrum 
analyzer. 
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In early phase tests the relative phase between unfiltered signals from the 
reference microphone and a selected microphone was measured using a Dranetz 
305- PA- 3009 phase meter. However, phase errors induced by the presence of 
harmonics in the signals complicated the measurement considerably. Consequently, 
later tests were performed using a Spectral Dynamics SD375 2-channel digital FFT 
spectrum analyzer which measured the transfer function between the two signals for 
a multitude of narrow passbands. Phase shifts of the fundamental component 
measured with the SD 375 are not subject to intrusion of harmonics and are 
therefore much more trustworthy. 

The control elements exist so that the experimenter need not continuously 
monitor the source section of the experiment. A signal is tapped from the reference 
microphone channel and is returned to the control input of the compressor. The 
compressor senses the rms voltage of the signal at the control input, compares it 
with a voltage set by the operator, and automatically adjusts the VGA gain so that 
the reference microphone signal stays constant. The compressor therefore 
guarantees a stable SPL at the reference position. 

C. Experimental Procedure 

The following is a description of a single segment of the experimental 
regimen as applied to an individual porous sample. A flowchart of the process is 
provided in Fig 4-4. 

1. Preparation 

Preparation for the experiments consumed nearly as much time as the 
experiment itself because the procedure of packing the tube was extremely tedious. 
We began by selecting one of three sample porosities, namely Q =0.94, 0.96, or 0.98. 
A large number of 3/4 in. squares were cut out of sheets of batted Kevlar^29. The 
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squares were originally in. thick and had a natural porosity of about 0.985. The 
squares were inserted one at a time a predetermined distance into the tube to 
create a sample having the desired porosity as defined by 

Kev 

where m is the mass of the sample, L is the total sample length, is the density 
of solid Kevlar (1450 kg/m^), and A is the cross-sectional area of the tube. The 
natural spring of the Kevlar resists packing, but friction with the tube walls holds it 
in place. Once the packing was completed the microphone holders and driver 
adapi“r were affixed to the tube. 

The homogeneity of the sample was tested by measuring the small-signal 
attenuation in the material at each of the microphone positions. With the source 
adjusted for a 100 dB SPL tone at 1 kHz at the reference microphone, the SPL at 
each of the microphone positions was measured. If the SPL differences between the 
signals at the 2, 5, 8 and 11 cm positions covered a range of more than 1 dB, the 
sample was pronounced inhomogeneous and rejected. In this case the hapless 
experimenter discarded the sample and repacked the tube. Otherwise, the propaga- 
tion tests could begin. 

2. Propagation Experiment 

The sound pressure levels of the signal harmonics at each of the six 
microphone positions were measured for all possible combinations of the following 
frequencies and reference microphone fundamental SPLs: 500, 1000, and 1500 Hz, 
and 100, 120, 140, 150, and 160 dB. It was not necessary to perform all 
combinations of microphone position, frequency, and SPL for the phase test. 
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Microphone selection was limited to the reference and the 8 cm and 11cm 
microphones because the changes in phase at high intensity were cumulative and 
were not fully developed at earlier microphone positions; at the 20 cm position the 
signal-to-noise ratio was usually too low to obtain stable phase readings. A full 
complement of reference SPLs was used, however, since the nonlinearly induced 
phase shifts tended to appear rather suddenly. Once the phase test was completed 
the microphones, microphone holders, and driver adapter were removed and the tube 
was reconfigured for the dc test. 

The accuracy of our measurements was potentially quite good. The B&K 
2010 is highly accurate and can be read to within ±0.2 dB. Digital sampling errors in 
the Spectral Dynamics analyzer give errors of ±0.2 dB and about 0.5° in phase. The 
HP spectrum analyzer, which was not designed for such exacting tasks, can be read 
to an accuracy of ±1 dB. The high level of accuracy is in reality unnecessary: 
inhomogeneities in the sample cause readings to deviate from simple small-signal 
attenuation by up to ±.5 dB Since the theoretical predictions assume a homo- 
geneous material, the inhomogeneitics are the limiting factor for the applicability 
of the data. 

3. dc Test 

The dc test was performed last for one reason: the sample usually 
compresses during the test, and the final porosity was unpredictable. Moreover, the 
sample can compress so far that it ceases to function as a termination. In extreme 
cases, the sample may recede beyond the last microphone. 

^Samples with larger deviations are rejected. 
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The test itself was quite simple. The flow rate was increased in steps 
that would be roughly equidistant on a log u plot. A large number of measurements 
were clustered at low flow velocities to ensure a good estimate of rj. As the flow 
rate was increased, care was taken to stay within the operating ranges of the 
flowmeters and pressure gages. The adapter block was occasionally removed for 
measurement of the sample length so that the compression effects could be 
accounted for at a later time. The experiment was completed when the range of the 
available pressure gages and/or flowmeters was exceeded. The sample was then 
discarded and the tube repacked. 

4. Compensation for Porosity Changes 

The purpose of the static flow resistivity test is to determine the 
resistivity of a sample as a function of flow speed. Implicit in this statement is the 
understanding that the frame is rigid, i.e., n= constant. Unfortunately, because the 
static flow test operates far below the decoupling frequency, the material is 
compressed by the flow, i.e., the sample length L is dependent on the velocity as 
L(u). As the sample compresses, the porosity is reduced and the resistivity of the 
material changes. These compressions do not take place for signals which, as in our 
experiments, lie well above the decoupling frequency. The subject of this section is 
a method of compensating the resistivity measurement for sample compression. 

The velocity dependent porosity S2(u) can be determined via Eq. (4.1) if the 

sample length L as a function of velocity is known. The sample length was measured 

for several velocities during the experiment, but the limited number of discrete data 

points obtained is hardly sufficient to define a function L(u). We therefore use the 

26 

method of cubic T-splines as developed by Foreman to generate an approximate 
curve from the discrete length data. The advantage of Foreman's method is that it 
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minimizes spurious curve "oscillations" and other inaccuracies caused by the 
inadequacy of a cubic polynomial fit. 

In order to distinguish between resistivity changes due to sample com- 
pression and those due to the expected amplitude dependence, we develop a model 
which describes the changes in resistivity with velocity and porosity. Hersh and 
Walker have determined the functional dependence of the linear resistivity rj on 
porosity to be 


r 


1 


(n) = -^ 
d g 



(4.2) 


where n is the dynamic viscosity, d is the fiber diameter, and g is an empirically 

determined constant of the material. From our data we calculate the average value 

1 

of g for batted Kevlar 29 to be roughly 0.065. Kuntz's data imply a value of about 
0.062, and Hersh and Walker reported a value of 0.059. We have used a 
compromise value of g=0.063 throughout these calculations. The nonlinear resis- 
tivity also changes with porosity, but its functional dependence has yet to be 
determined. We have observed, however, in Kuntz's data and our own, that the 
value of r 2 seems to be related to rj by a simple constant for a given material for a 
number of porosities. If we call this constant h, the resistivity can be determined 
from known values of rj,S2 , and h: 

r(fi,u) = rj(n) { 1 + h u ) . 

A more general form, which f^' ors neither Eq. (2.9) nor Eq. (2.10), is 

r(n,u) = r j(S2) f(u) , (4.3) 
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where the coefficients of the polynomial f(u) are porosity independent constants of 
the medium. For our experiment the porosity was a function of velocity 12 (u), so 
that our measured resistivity data is assumed to obey the relation 

r (Q(u),u) = r ,(12(u)) f(u) . 
rneas. I 

If the sample had remained rigid, the data would have obeyed the relation 

r(12(0),u) = rj(12(0)) f(u) , (4.5) 

where 12 (0) is the nominal porosity before any flow is applied. In order to adapt the 
data from a compressed sample to aft ideally rigid one, the two relations are 
combined to eliminate f(u). The rigid sample data is therefore related to the 
measured data by 


r 



rj (l2(0)) 


r 

meas. 



(4.6) 


Note that the explicit form of f(u) is unimportant here; it is sufficient merely to 
acknowledge its existence. 

D. Experimental Results and Comparison with Data 

This section is the confluence of the theoretical work of Chapter 3 and 
the experimental work of Chapter 4. Here we examine the correspondence between 
theoretical predictions and the experimental observations. The dc resistivity data is 
discussed first, since the resistivity properties of the medium are needed as inputs 
for the predictions. The acoustic data is then discussed in the same order followed 
in the theoretical discussions of Chapter 3: propagation of small-signal waves, 

finite-amplitude tones, and finite-amplitude distorted sinusoids. 


1. Resistivity 


Resistivity data for six samples, adjusted for sample compression, were 

2 

fitted to the u sgn(u) model, Eq. (2.9), and the u model, Eq. (2.10), by the method of 
curvilinear regression. The quality of fit, as measured by the mean squared 
deviation of the data from the fitted curve consistently favored the u sgn(u) model 
(see T ible B- 1 in Appendix B). A plot of some typical resistivity versus f 'ow 
velocity data is shown by the circles in Fig. 4-5 (see also Fig. 2-1). The data are 
somewhat scattered at low flow velocities because our equipment is less accurate 
for very small values of u and Ap. The calculated values of the constants rj and 
for the samples are shown in Table 4-1, along with the calculated values of the 
empirical constant g and the "relative nonlinearity" ratio The sample 

numbers were assigned in the order in which they were tested. 


TABLE 4-1 

Resistivity properties of porous samples 


Sample 

(nominal) 

■■l 

(mks rayl/m) 

■■2 2 
(rayl sec/m ) 

g 

rjCj/r 

1 

0.94 

50233 

6715 

0.066 

38.9 

2 

0.98 

10437 

3261 

0.071 

90.8 

3 

0.96 

29789 

6831 

0.062 

66.7 

4 

0.94 

51539 

12708 

0.064 

71.7 

5 

0.96 

28601 

6904 

0.066 

70.2 

6 

0.94 

52891 

11354 

0.062 

62.4 

7 

0.96 

28784 

6204 

0.065 

62.7 
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Data corrected for sample compress 
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illustration of the power of the compression compensation algorithm 
can be found in Fig. 4-5. The triangular symbols represent the raw data prior to 
correction for sample compression. The porosity of the sample remained constant 
at the value n=0.96 for flow rates below I m/sec. For higher flow rates, however, 
the sample compressed quite a bit; a flow rate of 1.73 m/sec reduced the porosity 
of this particular sample to f2=0.946. The uncompensated resistivity datum at this 
flow velocity was actually 62,139 rayl/m before the algorithm reduced it by 35% to 
40,538 rayl/m, where it falls right into line with the pattern set by the lesser flow 
rates. The remarkable consistency of the corrected points with those for which the 
sample was in fact rigid implies that the assumption underlying Eq. (4.3) i' valid. 
Adjustments made to resistivity data were similarly successful for all samples for 
which there was adequate L(u) information, save one. The exception was sample 2 
which, because of its very high porosity (12=0.98), compresse-' to about one third of 
its initial length. The cubic polynomial curve fit was simply insufficient to match 
the length data. 

Five of the seven samples studied had very consistent properties. Sample 
1, for no apparent reason, had a much lower relative nonlinearity than the others. 

2. Small-Signal Tests 

Attenuation and phase velocity predictions of Section 3.B were tested by 
measuring the change in SPL and phase of a propagating tone as a function of 
distance through the sample. The attenuation determined therefrom u plotted 
versus frequency and compared to predicted values in Fig. 4-6. The phase speed 
measurements and predictions are shown in Fig. 4- 7 . in both figures the solid lines 
represent predictions based on Eq. (3.23). The dashed line in Fig. 4-7 represents the 
isothermal sound speed value. Low frequency predictions agree well, but as the 
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frequency is increased the predicted phase velocity and attenuation curves diverge 
from data. Preliminary calculations indicate thnt most of the discre, .ncy could be 
removed by introducing a frequency dependent sound speed and a structure factor 
(see below). We have no explanation for the slight discontinuity in the data near 
1750 Hz in both figures. 

Recently published work by Lambert and Tesar^ apparently makes it 
possible to obtain an a priori value of the structure factor and the frequency 
ependent bulk modulus from resistivity tests no different from those discussed 
here. Lambert and Tescr's ti.eory is somewhat involved, but predictions of the 
small-signal phase speeds and attenuation constants for Kevlar'^'29 are accurate over 
a wide range of porosities, indeed, if *^ir work had been published earlier, their 
formulations of the structure factor and frequency dependent sound speed would 
probably have been incorporated in the present study. The behavior of these 
properties at high intensity remains, however, unknown. Nevertheless, predictions 
based on our simple theory are, by demonstration, reasonably accurate and have the 
advantage of being substantially less complex. A more detailed comparison of the 
two approaches, in:luding application to nonlinear effects, would be a good idea for 
a future study. 

The attenuation data serves as a check of the assumption that rj is the 
same for ac and dc flows. Table 4-2 compares the linear resistivity coefficients for 
ac and dc flows fo*" six of the samples tested. Data taken at low frequency (where 
our theory is highly accurate) yield estimates of the value of rj which are 
consistently lower than their dc counterparts by a few percent. The agreement is 
good enough, however, to validate our assumption that the ac and dc resistivities are 
equivalent for small signals. 


TABLE 4-2 


Comparison of measured static flow and acoustic resistivity 


Sample 

f2 

(nominal) 

r jdc 

(mks rayl/m) 

r jac 

(mks rayl/m) 

A% 

1 

0.94 

50233 

48831 

00 

• 

2 

C.^8 

10437 

9895 

-5.2 

3 

0.96 

29789 

29054 

-2.5 

4 

0.94 

51539 

48395 

-6.1 

5 

0.96 

28601 

27059 

-5.4 

6 

0.94 

52891 

50898 

-3.8 


3. Finite-Amplitude Tones 

We now discuss the SPL and relative phase data for finite-amplitude 
tones. The finite-amplitude signal at the reference microphone is not a pure tone, 
but rather a distorted sinusoid. Here we are concerned mainly with the fundamental 
and assume that the harmonics are small and have a negligible effect on the SPL of 
the fundamental. Theoretical predictions from Chapter 3 lead us to expect excess 
attenuation and a red'jction of the phase velocity at high amplitude. In other words, 
a finite-amplit' de tone behaves as a primary, as defined in Chapter 3. The 
nonlinear behavior is thought to be caused by the local increase of resistivity due to 
the passage of the wave. A good exam’ e of data that show excess attenuation is 
found in Fig- 4-8. The data is taken for sample 2, for which the porosity is 12=0.98. 
The three data sets correspond to three initial source amplitudes, and the dashed 
lines represent small-signal attenuation. No prediction curves have been included 
here because the tabulated value of V 2 for this sample is considered inaccurate. 
This data set contains the most pronounced excess attenuation we measured. The 
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Figure 4-8 
Excess attenuation 




explanation is that the lower small-signal attenuation rate allows the wave to 
remain strong for a longer period. The nonlinear effects are therefore relatively 
more important than in other samples. 

The other samples we tested showed very little excess attenuation (see, 
for example, Fig. 4-9) because their small-signal attenuation rates were large and 
our maximum source level was only 160 dB. (By using source levels of up +o 170 dB, 
Kuntz* was able to obtain large excess attenuations over a w’de range of porosities.) 
The three data sets in Fig. 4-9 once again correspond to three initial source 
amplitudes, but this time predicted propagation curves are included for each. At 
low intensity the predictions match the data well, but as the intensity is increased 
the expected amount of excess attenuation is exaggerated by the approximate 
solution. This is an example of the tendency of the approximate solution to 
"overshoot" as discussed in Chapter 3. Calculations have shown that the second 
derivative term discarded from Eq. (3.27) is moie important than had originally been 
thought. Because of the inaccuracy of the solution, we are unable to 
determine whether r 2 is the same for ac and dc flows. If the solution were exact, 
the validity (or lack thereof) of the measured static flow value would be immediate- 
ly apparent. Unfortunately, the solution is inexact. It was hoped that the 
measurements could be used as a test of the theory and, in particular, a check on 
the assumed value of r 2 , but since the solution of Eq. (3.29) misrepresents the 
theory, no direct assessment of the accuracy of the static flow of r,^ can be made. 
Nevertheless, the excess attenuation behavior is at least qualitatively predicted and 
is compatible with the concept of nonlinear ly augmented local resistivity. 

Measurements of the relative phase between two points in the wavefield 
indicate quite strongly that the phase speed of the pressure wave is a function of the 
wave amplitude. The relative phase is defined as the difference in phase between 
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the signals at tl»e reference and a downstream point, i.e., o A =®(x )-o(0). An 
increase of the relative phase corresponds to a decrease in the phase velocity. 
Figure 4- 10 is an example of some typical data. The solid curve represents 
predicted values of the relative phase as a function of SPL. Abscissa values 
represent deviations of the relative phase from t'le predicted small-signal value 
A0ss» that is, 

'^»Sraph= • 

Vertical offset errors due to inaccurate predictions of qj^ are present which 
complicate the picture somewhat. Hov'ever, the shape of the curve is a generally 
dependable measure of the change in phase velocity. 

The measured relative phase data did some things we expected, and others 
we did not. As the SPL was increased, the relative phase usually decreased at first, 
but always increased at higher levels. The decrease took many forms; sudden jumps 
were as common as gradual changes. The overall tendency of the relative phase to 
increase with SPL is compatible with the concept of augmented resistivity, because 
the data show that the phase velocity of the pressure wave was reduced at high 
amplitude. The decrease in the relative phase was totally unexpected. Small 
hysteresis loops were also observed; that is, the data followed a different path 
depending on whether the SPL was increasing or decreasing. The anomalous 
decrease, sudden jumps, and hysteresis behavior suggest a nonlinear bistability for 
which the local stocfiastic variation of the quantities tj and r 2 causes there to be 
more than one possible wave number at high intensity. 

The skeptical reader might suspect that the hysteresis behavior was 
introduced somewhere in the measurement apparatus. The measurements were 
repeatable, however, for various microphones, microphone positions, frequencies. 
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and even different FFT analyzers. We therefore believe the nrieasure'nents to be 
true representations of actual acoustical events. 

We conclude from the finite -amplitude tone measurements that atten- 
uation is indeed increased and phase velocity reduced at high amplitude. The 
approximate solution gives good qualitative predictions for both. Although we 
should have liked to draw a conclusion concerning the values of r 2 for ac and dc 
flows, the inherent inaccuracy of the approximate solution prevents us from doing 
so. We are ther^iore unable to confirm or deny our initial assumption that the dc 
and ac resistivity are interchangeable at all particle velocities. 

4. Distorted Sinusoids 

Having finished with the fundamental, let us now concentrate on the 
higher harmonics of the signal. The higher harmonics generated by the propagating 
fundamental are no longer neglected, nor are those present at the reference 
position. Analyzing the behavior of the higher harmonics gives us information about 
the form of the r(u) function. A computer program was written to numerically solve 
Eq. (3.29) for periodic waves. A cubic distortion pattern is predicted and the 
harmonic propagation! curves are predicted to cross. Figures 4-11, 4-12, and 4-13 
give measured and predicted harmonic propagation curves for three different 
samples. In each case, data for the fundamental show that it behaves essentially as 
treated in the previous section. In fact, the small difference between the analytical 
(see Fig. 4-9) and numerical solutions for the fundamental demonstrates the minimal 
effect of the higher harmonics on the fundamental. The SPL of the second harmonic 
is seen to decrease monotonically with distance just as the SPL of the fundamental 
does. Therefore the second harmonic also functions as a primary. The third 
harmonic SPL, however, increa*-,es initially and routinely equals or exceeds he level 





of the second harmonic. The upshifted energy from the fundamental skips over the 
second harmonic and reinforces the third. As a result, the propagation curves for 
the second and third harmonics may c ross. This is a distinguishing mark of the cubic 
distortion pattern which was ('iscussed in Chapter 3. An ordinary quadratic 
nonlinearity, such as that for hydrodyna 'c distortion or that assumed by Kuntz, 
would cause the second harmonic to beha/e as a distortion product rather than a 
primary. A quadratic nonlinearity does not therefore fit the distortion pattern 
observed here. Data for the fourth and fifth ha'^monics indicate that both receive 
significant energy from the primaries. Their behavior, however, neither confirms 
nor denies a cubic distortion mechanism. 

The predicted harmonic propagation curves art in substantial agreement 
with the data. Predictions of the SPL of the fundamental and second harmonic are 
too low because the numerical solution exaggerates the nonlinear effects. The 
overshoot phenomenon has a different effect on distortion product predictions. 
Since the solution exaggerates the nonlinear effects, the approximate solution icr.Js 
to overestimate the amount of energy received by a distortion product and thus 
overestimates the SPL as well. The error is balanced somewhat, however, because 
the approximate solution also underestimates the strength of the primaries which 
drive the distortion products. 

Prediction curves for the higher harmonics, especially those that receive 
a great deal of upshifted energy, often have cusps (sc*e n=5 curve in Fig. <f-ll). The 
cusps are artifacts of the solution procedure. The cause lies in the fact that the 
differential equation is solved for initial conditions at the reference postion X =0. 
The inhomogeneous part of the solution, representing energy upshifted from the 
primaries, is set to zero by the computer program at this point. In reality, the 
upshifted energy is nonzero at the reference position: it is only zero at the true 
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physical boundary, the material/air interface, which is not located at X=0. The nth 
component at X=0 is partly due to energy upshifted while traveling from the physical 
boundary ax the reference position. Since the waves represented by the two parts of 
the solution travel at different speeds^ the upshifted signal is occasionally out of 
phase with the signal represented by the homogeneous solution. Destructive 
interference takes place, and one and occasionally more spurious cusps may appear 
in the predicted propagation curve for a harmonic. 

E. Summary of Experimental Work 

As a test of the validity of our theory, we have performed a series of 
experiments to examine the behavior of acoustic signals in batted Kevlar 29 at low 
and high intensity. The following five conclusive results have been obtained . 

1) The static flow resistivity depends on the flow speed as 
r=rj+r2usgn(u). 

2) The value of rj for oscillatory flows is equivalent, or nearly 
equivalent, to that for static flows. 

3) A finite amplitude tone undergoes excess attenuation which is 
somewhat less than predicted. 

4) A finite-amplitude tone propagates more slowly than its small signal 
counterpart. 

5) The harmonic distortion pattern has cubic character. 

2 

Although point (5) admits both the u and sgn(u) resistivity models, the evidence 
from (1) implies that the sgn(u) model is correct. 
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Consider that the harmonic tends to iiravel at C 


■pn 


comes the inhomogenous part, trave' at Cpj and Cp 2 * 


while the primaries, whence 


In addition, two results have been obtained which, although not con- 
clusive, are considered trustworthy and merit further study: 

1) the resistivity can be separated into functions of porosity and flow 
velocity as r=rj(fi)f(u), where the as yet unspecified coefficients in 
the function f(u) are independent of porosity, 

2) an apparent nonlinear bistability has been observed which indicates 
that the local stochastic variation of the material properties may be 
of importance at high intensity. 

We have not been able to deternine whether the value of C 2 is the same 
for ac and dc flows because of the inaccuracy of our solution. 
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CHAPTER 5 


CONCLUSIONS FOR CHAPTERS 1 THROUGH 4 

The propagation of finite amplitude sound in rigid air-filled porous 
materials has been studied both theoretically and experimentally. In particular we 
have studied the effects of amplitude dependent resistivity on propagation of 
acoustic signals. Expeiirnents performed on batted I'evlar 29 are in at least 
qualitative, and in many cases quantitative, agreement with the theoretical predic- 
tions. 

The theory is based on the following model of the pKjrouo material: the 
material is rigid, incompressible, and homogeneous and has only two important 
properties, porosity fl and resistivity r. Detailed knowledge of the material 
geometry is unnecessary for this model. The resistivity has been measured for dc 
flows and is shown to be the following function of flow velocity u: r=r j+r 2 usgn(u) 
(see Eq. (2.9)).^ Previous work by other authors shows that the static flow 
resistivity function may Le applied to oscillatory (acoustic) flows. The components 
of the flow resistivity play two different roles. The "linear resistivity" rj 
determines, among other things, the small-signal attenuation, which in a typical 
porous material is quite severe. The amplitude dependent resistivity leads to strong 
nonlinear effects, namely excess attenuation and harmonic distortion. The effects 
of hydrodynamic distortion, which normally give rise to shock formation, are 

t 2 

The model of Eq. (2.10), r=r j+r^u , is far from useless; a discussion of results based 
on this model is given in Appendix B. 
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completely overshadowed by both the linear and nonlinear resistivity effects and are 
neglected. The isothermal sound speed is used ior all frequencies. We have 
therefore established the resistivity as the prime mechanism for attenuation and 
distortion of acoustic signals in porous materials. The apparent lack of sophis- 
tication of the model respresents a deliberate retreat from more elaborate theories; 
the simplified theo'^y has the advantage of being easily extended to include finite 
amplitude behavior. 

A nonlinear wave equation, Eq. (2.23), is derived from the model. The 
wave equation has a relatively uncomplicated form, but turns out to be very 
difficult to solve in the general case. Transformation to the frequency domain 
produces an infinite set of coupled inhomogeneous Helmholtz equations, one for 
each harmonic. 

In the low intensity limit the Helmholtz equations are homogeneous and 
solutions are easily obtained. Because the propagating acoustic wave is subject to 
diffusion (see Chapter 2), we expect to observe high attenuation rates, severe 
dispersion, and low phase speeds. 

Solutions for high intensity waves are considerably more difficult to find. 
An approximate first integral, Eq. (3.29), of the Helmholtz equation was derived 
because numerical calculations for the solutions of the second order Helmholtz 
equations (Eq. (2.26)) were unstable. Two different approaches were used to solve 
Eq. (3.29). In one the higher harmonics were neglected entirely and in the other the 
higher harmonics are taken into account. Using the first approach we are able to 
obtain an analytical solution from which expressions for impedance, and the 
amplitude and phase speed of both the particle velocity and pressure waves are 
found. Ml of the many informative results obtained from these solutions may be 
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summarized in the statement that the mean resistivity increases with the wave 
amplitude. Finite-amplitude tones are therefore predicted to have higher atten- 
uation rates and impedances and lower phase speeds than small-signal waves of the 
same frequency. Saturation behavior is predicted at very high intensities; i.e., 
nonlinear losses determine the maximum pressure amplitude which can be reached 
at a given point in the wavefield. 

Numerical solutions of Eq. (3.29), including all harmonics, are performed 
by a computer program written for this purpose. Predictions based on this solution 
lead us to expect a cubic harmonic distortion pattern in which the fundamental 
energizes the third and higher odd harmonics. The second harmonic is predicted to 
behave as a fundamental as well and to directly energize every other even harmonic, 
that is, the sixth, tenth, fourteenth, etc. The fourth harmonic and multiples thereof 
are products of the interaction between the first and second harmonics. 

Several experiments were performed to investigate the propagation of 
finite-amplitude sound in porous materials. A traveling wave tube filled with a 
porous sample was fitted with six microphones, a compression driver, and extensive 
electronic devices for the acoustic measurements. The tube, with the same sample 
intact, was reconfigured with pressure gages and flow meters for static flow 
resistivity measurements. The samples were batted K.evlar^29 and had porosities 
fl=0.94, 0.96, and 0.98. Smell-signal measurements were taken in the frequency 
range 400-6200 Hz. Finite amplitude measurements were conducted at a frequency 
of 1 kHz. Sound pressure levels of up to 160 dB (at the reference position) were 
used. 

Experimental results were in general consistent with predictions. 
Measurements at low intensity showed the attenuation and phase s[>eed were strong 


107 


functions of frequency. Measured attenuations for all samples were of the order of 
hundreds of dB/m and increased from low values at low frequency towards 
asymptotic values at high frequency. Phase speeds were substantially reduced at 
low frequency and increased to near the isothermal sound speed at high frequency. 
At high intensity, measurements of the SPL of the fundamental with distance 
showed a small but significant amount of excess attenuation. Because of the limited 
power-handling capacity of our acoustic source, however, the saturation region was 
out of reach. The approach to saturation has been much better documented by 
Kuntz. ^ The relative phase of the fundamental between the signals at two positions 
was measured as a function of amplitude. An increase in relative phase corresponds 
to a decrease in the phase speed. The measured relative phase was stable at low 
intensity, usually decreased slightly as the amplitude was increased, and always 
finished with a strong increase at the highest SPLs tested. The decrease was not 
always gradual, but often involved sudden jumps of several degrees for a source 
level change of 1 dB. In addition, small hysteresis loops were observed in the 
relative phase data. Finally, the harmonic distortion pattern was indeed observed to 
be cubic, that is, the first and second harmonic both behaved like "fundamentals"; 
that is, they decreased monotonically in strength with distance. The third harmonic 
SPL, by contrast, was observed to routinely exceed the amplitude of the second. 
The fourth and fifth harmonics often increased initially in strength as well. Such a 
distortion pattern could not be caused by a hydrodynamic distortion mechanism and 
is solely attributed to nonlinearity of the porous material resistance. 

We now compare the theoretical predictions and experimental obser- 
vations with a view towards 'Evaluating xhe validity of the model. Excellent 
qualitative agreement was obtained tor all predicted pheromena for wh»ch 
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measurements were attempted. Quantitative agreement ranged from satisfactory 
to excellent. 

(1) Low intensity attenuation and phase speed measurements were in 
very good quantitative agreement with predictions based on the 
measured dc value of rj; at low frequency (< 1500 Hz) the agreement 
was excellent. We therefore surmise that the dc and ac values of rj 
are the same. 

(2) Measurements of sound pressure level versus distance at high 
intensity showed some excess attenuation, but in general we 
observed much less than was predicted. The assumptions on which 
Eq. (3.29) is based were apparently more easily violated than we had 
expected. It was hoped that this solution could be used to 
demonstrate the equivalence of r 2 for ac and dc flows, but the 
inaccuracies prevent us from drawing a solid conclusion. Never- 
theless, the dc value of r 2 seems to be at least approximately 
correct for use with acoustic signals. The relative phase data 
deviates from predictions as a result of a process which we suspect 
is a nonlinear bistability induced by the inhomogeneity of the 
material. The predicted increase in relative phase (a decrease in 
phase speed) appears, however, consistently at high intensity. 

2 

(3) The existence of cubic distortion is compatible with both a u sgn(u) 
or u^ nonlinearity. The gross details of the harmonic propagation 
curves are very adequately predicted; for the most part the pre- 
dictions are within 5 dB of the measured data. 


On the basis of the qualitative and quantitative agreement between the 
theory and experiment, we conclude that 

1) the resistivity is the dominant mechanism for acoustic loss and 
distortion in a porous material, 

2) the mathematical form of the resistivity function is r=r j+r 2 usgn(u) 
for static flows and is retained for acoustic flows, and 

3) the value of rj is indeed the same for ac and dc flows, but 

4) the agreement between the value of r 2 for ac and dc flows cannot 
be determined because solutions are inaccurate for finite-amplitude 
waves. 

Two phenomena were observed which suggest modifications for future 
theories. Small-signal measurements indicate that the accuracy of the theory would 
benefit somewhat from a more sophisticated model in which the sound speed were 
allowed to be frequency dependent, and a structure factor were included. However, 
there is currently little knowledge available on the variation of these parameters at 
high intensity. The second phenomenon is the hysteresis loops that were observed in 
the phase speed versus amplitude data. The nonlinear bistability implies that local 
stochastic variations of the material resistivity, i.e., material inhomogeneity, may 
be important at high intensity. 

The experiments had two basic shortcomings First, the static flow test 
compressed some high porosity samples so much that even our otherwise highly 
effective compression compensation algorithm could not account for the changes. 
Use of lower porosity samples solved the first problem but led to the second; the 
source amplitude was not high enough to observe significant excess attenuation for 
all the porosities tested. Results of future experiments will be more conclusive if 
these two problems can be solved. 
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Three topics are suggested for future study. First, there is a definite 
need for more accurate solutions of the infinite set of coupled second order 
inhomogeneous ordinary differential equations represented by Eq. (2.26). Second, 
the accuracy of theoretical predictions may benefit somewhat from a more 
elaborate model of the material which includes structure factor and a frequency 
dependent sound speed. Furthermore, the change in these parameters at high 
intensity should be studied. Third, since the material is likely to be used in an 
inhomogeneous state, the effects of the inhomogeneity should be investigated. 

The theory developed here was originally intended to be applied to the 
problem of absorption of high intensity sound in a lined duct. Theoretical work on 
this topic is already being pursued by this author. I. point- and local-reaction 
models of duct liners are to be used with confidence, it will be necessary to 
establish the validity of amplitude dependent impedance expressions derived in this 
work. Despite the difficulties encountered, the theory we have developed appears 
to offer the best hope for solving the important and practical problem of high 
intensity sound absorption in lined ducts. 
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CHAPTER 6 


THEORY OF LINED DUCTS 

A. Introduction 

Chapter 6 is an application of the accumulated knowledge of the fore- 
going chapters to the problem of a lined duct. We are particularly interested in the 
effect of the materia! nonlinearity on the absorptive power of the duct. We begin 
by deriving a nonlinear nondimensional vector wave equation for acoustic waves in 
the material. Next we investigate the process of nonlinear reflection from a porous 
half-space. Lastly, a third order perturbation solution for a lined duct is prepared, 
but is never executed because of its forbidding complexity. Instead, we propose an 
ad hoc model to aid in understanding the nonlinear effects. 

B. Nonlinear Nondimensional Vector Wave Equation 

In previous chapters the nonlinear resistivity was calculated assuming 

one-dimensional flow. In order to consider a lined duct, however, a general 

expression for the resistivity is needed. In particular, we are interested in finding 

vector expressions to replace Eqs. (2.9) and (2.10) for use in a vector wave equation. 

2 

Of the two previously used forms ci the nonlinear resistivity u sgn(u) and u , th;. 
former represents the instantaneous speed of the flow aind the latter the square of 
the velocity. It turns out that the square of the instantaneous particle velocity can be 
expressed in vector terms as u»u =|u|^. The vector form of Eq. (2.10) therefore 
lends itself rather well to a perturbation method of solution. The instantaneous 
particle speed is of course >/u*u . An expression of the resistivity based on this form, 
however, does not allow a perturbation solution because the square 
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root function is not distributive. Even though we prefer the sgn(u) formulation (as 
discussed extensively in previous chapters), we use the quadratic form of the 
resistivity in preparation for a perturbation solution. 

Since the porous material is assumed isotropic the one-dimensional 
particle velocity u can be replaced by its vector counterpart, 

p^ > Pq y • u = 0 . (b.l ) 

The nondimensional continuity equation, expressed in vector notation, is 

p^ + V . V = 0 , (6.2) 

and the momentum equation becorr.es 

+ RjV + V • V )) =0 , (6.3) 

where the vector form of the quadratic nonlinear resistivity R^v*v has been used. 
The wave equation is found by adding the results of the operations - y/Q on the 
continuity equation and dfd^ on the momentum equation; 

v^T - — > RjVt » R 3 (v(v = 0 . (6.4) 

This equation is awkward in vector form and has the added disadvantage of allowing 

rotational solutions. Since we are interested in scalar fields describing irrotational 

motions, two modifications are made: the velocity vector is first replaced with the 

gradient of the scalar potential field V 0 , 

2 

(v0(^0 • ^0)^ = 0 (6.5) 

and the divergence of the result is taken to reach the scalar form 

{0^^. V ^0 + Rj 0 ^) + Rj^(V ^0 +V 0 . V )(^0 . ^ 0 )j = 0 . 


( 6 . 6 ) 



The frequency domain version of this equation is found by substituting for 


the potential in the by now familiar manner to give 


vW ♦ = -jnn R3EE p* ^ Vp-q^ ’ 

p q ^ ^ ^ 


where 


and 


(X) exp(jnt) 


The perturbation solution is prepared by ex, es'-In,, the nth frequency 
component of the potential field as 


0p, = ^ ^n3 (6.8) 

where only the first and third orders are represented because the second is 
I'lomogeneous and therefore redundant with the first. The value of 4 will be 
specified later. Substitution into Eq. (6.7) and grouping b> powers of epsilon gives 
the following relations: 


V2(r^qn^>«’nl 


(6.9) 


and 

vV > qn^)0p3 =-jnnR3SL<^\l " * ^^n-p-ql^ * 

Once the first order field is sp>ecified, the third order field can be found by 
operational integration. 

Later we find that the p>erturbation technique is very awkward for 
application to the problem of li led ducts because of the enormous complexity of the 
nonlinear interaction of the wave fields. The perturbation solution nevertheless 
serves to illuminate several interesting physical processes including 


r< i lotion of an intense wave from a porous half -space as discussed on the next 
section. 

C. Reflection from Porous Boundary 

1. Small-Signal Problem 

The problem n{ the reflection of sound from a porous material is 
interesting for low intensity sound waves, but the inclusion of the effects of 
material nonlinearity leads to some surprising phenomena. It turns out that at high 
intensity the wave reflected f -om a porous surface does not have the same spectrum 
as the corresponding incident wave. This prediction relates directly to experimental 
obst vations made by Zorumski and Parrot. The following investigation is 
predominantly theoretical in nature. 

The wave equation is very complicated, but still lends itself to a 
workable, if somewhat tedious, perturbation solution. To begin the problem, we 
first define two potential fields: 0^ in the air and 0^ in the material. They of 
course have the frequency domain manifestations 0^ and 0^. These fields are 
further broken down as 

0"^ = €0^, + ■>■ . . . (6.11a) 

n nl n3 

and 

0n = 

where € is the magnitude of the potential function of the incident wave. Since 
nonlinear effects in the air are neglected, the Helmholtz equation for the air is 
homogeneour. We shall see, however, that it is stilt necessary to retain the third 
order solution in the air because it represents ene. gy which is ladiated from the 
materi." surface. The first and third order Helmholtz equations are: 


and 


OnVo 

0 , 


A 

h3 


( 6 . 12 ) 


(6.13) 


-A ^A ^A 
®nk = ^nkl * ^nkR 


In the material, Eqs. (6.9) eind (6.10) apply. To further complicate matters, the wave 
fields in the air are composed of incident and r ‘fleeted waves, 

(6.14) 

The material is assumed to be of semi-infinite extent, so that the wave fields are 
composed '•olely of transmitted waves, 

C-<T • 

The conditions that pressure and mass flux must balcmce are invoked at 
the boundary. The pressure boundary condition is written in potential form from the 
equation of continuity: the substitution of V0 for 7 yields 

p +V^=0 (6.16) 


and in the frequency domain this expression becomes 

inP = 

' n ^n 

equating the pressures on either side of the boundary at X = 0, we find that 




X=0 




x=0 


(6.17) 


(6.11i) 


Since the density of the air does not change across the boundary, the boundary 
condition for the conservation of mass is simply that the volume velocities match. 
In the air the volume velocity is 

kA 


90 , 


nk 


ax 


X=0 
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and in the material 


« A 


ax 


X=0 


where A is the cross-se<'tional area of the boundar y. The second boundary condition 
is therefore 


9<^nk 

ax 


x= 0 


n 


ax 


X = 0 


(6.19) 


The first order wave fields are defined as follows: 


and 


{'Jy 


nil = \l H W ^ 


n<>) . 


«>nm = Bnl '‘P (-ijf * * »‘"*n « >) 

®MT = C„,«p(-)(q„xX *qn£«>) • 


( 6 . 20 ) 


( 6 . 21 ) 


( 6 . 22 ) 


2 2 2 

where = q^ . The situation is depicted in Fig. 6-1, where X is the coordi- 

nate normal to the surface and 4 is the coordinate parallel. The algebraic 
manipulations encountered in applying the boundary conditions are tedious and not 
very informative; hence we simply show the results. The small-signal pressure 

reflection coefficient pAo ) is equivalent to B ,/A ,, which is 

n n ^ nl nl 


p (a ) = 

*^n n 


^ ^ Jy% 


(6.23) 


The pressure trcinsmission coefficient t^(fl^) is simply equivalent to one plus the 
reflection coefficient or 
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Figure 6-1 

Reflection and transmission of 
small -signal wave having rich spectrum 


(6.?4) 


(e) = 


‘2'fy % cose^ 




nx 


The value of C ,/A , is not equal to tj.0 ), however, but 
ni ni n n 


j\L - rnL. t(e) 

A 7 n n 

vq. 


(6.25) 


The magnitude of the transmitted waves, whence comes the nonlinear part of the 
solution, is determined by ihe value of 

The angle of the phase fronts within the material is defined by 


ib^ = tan 


-1 


Re 


fte 

■^nX 


(6.26) 


where 


and 


q„X =Vqn 


(6.27) 

It is important to note that the angle of travel within the material is frequency- 
dependent for all angles of incidence except normal incidence. Since low frequency 
waves have lower phase speed in che material (see Fig. 4-6), they are refracted 
strongly. High frequency waves, on the other hand, are refracted less. An incoming 
signal with a rich spertrum therefore fans out into its respective components as a 
consequence of the severe dispersion in the material. The complete first order field 


in both the air and the material can now be defined from the values of A„, and 

ni n 
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2. Nonlinear Reflection and Refraction 


In this section we seek practical insights into the problem of reflection >f 
high-intensity sound from porous materials. This work is done in preparation for the 
solution of the lined duct problem, and to help us understand some of the vagaries of 
that problem. Here we merely prepare the theoretical tools necessary for explicit 
solution of reflection problems; we choose to concentrate on explaining qualitative 
results for instructive cases. 


We first define the third order wave fields for the 
transmitted waves: 

reflected and 

0 n3I = ° ’ 

(6.28) 

0n3R = ®n3 ®n3 ^ ^ ®n3^ ’ 

(6.29) 

and 

«’nJT= Cn3'’‘p(-)V3**'’n4 3«>) • 

(6.30) 

where, for a nonlinear interaction involving the rth, sth, and n-r 
components, the following relations apply: 

-sth harmonic 

*^nX 3" ^rx * %X * ^n-r-sx ’ 

(6.31) 

^n^ 3“^r4 ^ %-r-s4 * 

(6.32) 

®n3 ~ [o/n) (r sin 0^ + s sin + (n-r-s) sin 0^ ^ ^)j 

, (6.33) 
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and 




(6.34) 


The third order incident wave field is identically zero because the nonlinearity of 
the air has been neglected. The subscript 3 attached to variables in the following 
discussion indicates that they apply to the third order solution. In many cases we 
will study, the travel directions and wave numbers for the third order waves are 
redundant with those of the first; but in order to be sufficiently general, the explicit 
expressions are given. 

It is clear from Eq. (6.10) that the third order transmitted field is nonzero 

for at least some harmonic components of an intense wave. Since the same 

boundary conditions apply for all orders of the solution, it is necessary that a 

component generated in the material by nonlinear effects must be balanced by a 

component in the air, which then radiates from the surface. This means that the 

reflected wave field carries with it the imprint of the nonlinear interaction which 

27 

takes place inside the material. Observations made by Zorumski and Parrot 
corroborate this prediction; while studying the reflection and transmission of sound 
through thin porous sheets, they observed that the reflected wave field for an 
incident pure tone contained a strong third harmonic. 

The first nonlinear case we consider is that of an intense pure tone 
incident on the material at some angle 6j, as depicted in Fig. 6-2.^ The transmitted 
fundamental follows the angle which is determined from Eq. (6.26). The third 
order field contains two components; a correction to the fundamental proportional 


Circled numbers denote intense harmonic components. 
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lines of equal onplitude 



Figure 6-2 

Reflection and transmission 
of intense pure tone 


2 3 

to 0 j and a third harmonic distortion product which comes from 0j. Distortion 

products such as these will be referred to as "products of self-interaction," because 

the nonlinear interactions by which they are produced involve only one frequency 

component. It is interesting to note that these distortion products travel in the 

same direction as the primary (in this case the fundamental) from which they 

originated. It is easily shown that 

'^3X3" ^^1 ‘^3^3" ^^1 therefore = 'tj 

and 


^1 X 3 “ ^1 '^l 4 3 " ^1 therefore 'I'j ^ = 4'j . 


This leads to a second observation: products of self-interaction travel collinearly in 
the material with the primary from which they are generated. It is also of interest 
to note that the angle * 'I'^. In other words, the nonlinearly produced third 
harmonic travels in a different direction from one that originates outside the 
material. Furthermore, the angle which the distortion products take on emerging 
from the material is equivalent to that for the emerging primary which energized 
them. This is easily seen from Eq. (6.33). If the solution were carried to higher 
order, all of the odd harmonics would be represented in the transmitted and 
reflected waves, and all would be collinear within and without the material. 

The reflection coefficient of the fundamental is affected strongly by the 
material nonlinearity. In Eq. (3.41b) the predicted dependence of impedance 
magnitude on amplitude is given. Figure 6-3 shows the predictions of this equation 
and data obtained from Kuntz,^ where the value of K has been selected to best fit 
the data. Our predictions are qualitatively borne out. With this in mind, we use 
Eq. (3.40b) using the fitted value of K to determine the normal incidence energy 
absorption coefficient as a function of SPL. Figure 6-4 is an illustration of the 
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batted Kevlar 29 



Normal incidence energy 
absorption coefficient versus SPL 






severe dete-ioration of the absorber efficiency as the signal strength increases. 

The second nonlinear case we discuss is that of an incident wave made up 

of two intense components (primaries); the fundamental and the second harmonic. 

Once again, these primaries travel into the material in tne directions designated for 

linear behavior (see Fig. 6-5). This time, however, there are six components in the 

third order field. Each of the first six harmonic components i: represented. The 

first and third qualify as products of self-interaction, and travel collinearly with the 

fundamental. Likewise, the second and sixth harmonics are direct products of the 

second harmonic, and travel collinearly with the second. The fourth and fifth 

harmonics, as products of a mutual interaction, do not travel either of these two 

paths; they travel instead in intermediate directions. These directions can be 

2 

calculated rather easily. The fourth harmonic is a result of the product 0j02» and 
has the wave numbers 

^4X 3 = ^‘’iX ^ ^2X 
and 

which gives the propagation direction 


sin e 


♦ 43 = tan 




2 

Likewise, the fifth narmonic is a result of the product 02 0j, and has these 
wave numbers and propagation directions 


*^5X 3 " ^IX ^ ^'^2X 

'^54 3 = ^14 + 2q2^ =5/v/ysin0j 
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_j/5/Vy sinflA 

*53 = l^Rc[q5,3]-j • 

Perhaps the most interesting point to note here is that despite the unu'^ual 
direction of travel for the harmonics (relative to their linear paths) their reflected 
wave counterpa’ts are still collinear with the first order reflected wave. This is 
easily verified from Eq. (6.33). This brings us to another principle: if th“ incident 
signal is collinear in the air, then neither the frequency-dependent refraction nor 
the nonlinear interaction between the nonrollinear elements in the material 
prevents the reflected wave field from emerging in a collinear fa*i'lon. 

A similar example is given in Fig. 6 - 6 , where the case of an intense first 
and third harmonic incident on the material is shown. Qualitatively, this case is no 
different from the previous one, except that the third harmonic distortion product 
travels in a different direction from the incident third harmonic. The result is 
spatial interference patterns resembling Moire patterns. 

The final example we consider is that of an intense first and second 
harmonic incident on the material, but this time at different angles and $ 2 ’ Ry 
all appearances Fig. 6-7 is no different qualitatively from the case in Fig. 6-5, but 
an important distinction appears in the reflected wave field: some of the wave 

components are rxj longer collinear with the incident primaries in the reflected 
wave. The reason for this is not difficult to .hew. The {-direction wave numbers 
for the fourth and fifth harmonics are 

q«l 4 3 ' + sin 02 ) 

and 

^54 3 = (i/'/y) (sin^j + a sin 62 ) • 
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The emergence angles are therefore equal to 


and 


I / sin flj + sin 02 ^ 
e„=sin- I 5 


= sin 


-1 I 


+ 4 sin 02^ 


We see that if noncollinear incident components interact within the material, the 
products of their mutual interaction in the reflected field will not propagate in 
either of the primary directions, but will adopt instead some intermediate direction. 

The purpose of this discussion has been to understand in a practical way 
the nonlinear interaction in the material, and its effect on the wave fields inside and 
outside the material. In particular we have been interested in the wave numbers and 
propagation directions of tlie nonlinearly induced distortion products. The fact that 
signals can arrive from and depart in any direction from the material creates no 
difficulties in the reflection study. Difficulties arise, however, in the lined duct 
problem because wave numbers and propagation directions are rigidly linked by the 
waveguide. The nonlinearly induced waves in the waveguide guarantee that there is 
a certain amount of mode conversion. The resulting algebraic tedium of the 
perturbation solution of the lined duct problem renders it nearly useless. 


D. Lined Duct Problem 

Our original plans to solve the lined duct problem by perturbation were 
very ambitious. We first intended to solve the small-signal problem of a nonlocally 
reacting liner. The theory was then to be expanded to include the nonlinearity of 
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the duct liner, and lastly to include the nonlinearity of the air as well. Unfortu- 
nately, this plan had to be curtailed considerably in the middle of the first phase: it 
is difficult or impossible to obtain general solutions for the modal wave numbers. 
After making some approximations, we arrive at an ad hoc model for calculating 
attenuation of a nearly-plane wave having wavrlength much larger than the duct 
dimensions. A computational scheme is proposed to visualize the effects of liner 
nonlinearity on the absorptive efficiency of the duct. A computer algorithm is also 
suggested to calculate the optimal duct liner configuration. 

We first define the waveguide geometry. The reader will recognize in 

Fig. 6-8 the axes from the reflection problem. The following derivation conforms 

28 8 

roughly to those performed by Kurze and Ver or Scott. We define two 
wavefields, consisting of a wave traveling in each direction in both the air and 
material: 


and 


where 


0^^ =(g'^ exp(-ja X) + exp(ja xA exp(-ib 4) 

^ mn y mn ' mn mn mn J ^ mn 

0'|'_ = exp(-jc x) + exp(jc x)\ exp(-jd 4 ) , 

mn \ mn ’ mn mn ^ ’ mn j ’ mn 


(6.35) 


(6.36) 


The wave numbers and modal solution will be defined by the previous two boundary 
conditions, Eqs. (6.18) and (6.19), plus 


and 


30 ‘ 


mn 


ax 


= 0 


X=-L; 


(6.37) 
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= 0 . 


(6.38) 


The two new conditions represent the requirement that there be no normal velocity 
component at the walls. Upon applying the boundary conditions, we find the 
following equation for the wave numbers 


mn mn A 


where 


“"''mn'-vi* ’ 


(6.39) 


a - c = q 

mn mn y ^n 


(6.40) 


and b = d .A major complication arises at this point because of three factors; 
mn mn 

tangents occur on both sides of Eq. (6.39), and the wave numbers are both complex 

and hyperbolically related. To our knowledge, the problem of determining the wave 

numbers has not been solved satisfactorily because numerical solutions of these 

simultaneous equations are extremely unstable. Standard practice seems to be to 

assume something about the duct liner, for instance, that it is locally reacting (see 
29 

Cremer ). This simplifies the mathematics considerably. It becomes necessary at 
a later stage for us to make approximations as well, but for the moment we assume 
that the wave numbers can be found, and proceed with a third order solution for the 
duct. The following four wave fields are defined: 


®mn3H m3H *’‘P<-i“m3’' > * m3H '’‘P<i^m3 * > ) « > 


'”‘P'-i'>mn3 « > 

(6.41) 


(6.42) 
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®™3p = (° mSp *> * m3p '’‘P<i'm3 ’<>) '’‘P<-'‘‘mn3 « ’ 

®mn3H =('^“ 3H '’'P<-''mn3 * > * = m3H '’“P<i''mn3 *>) « > 


where the two homogeneous fields ( 0 fpri 3 }^) exist for the purpose of satisfying the 

boundary conditions at 4=0. The wave numbers of the waves radiated into the air 

space are not generally those which fit the modes of the guide, hence this set of 

homogeneous modes is needed to balance the nonlinearly induced ones. 

The harmonic interactions are calculated from Eq. (6.10), and expressions 

for the wave numbers and mode strengths are found by applying the boundary 

2 2 

conditions. The number of calculations is potentially large: there will be 9M N 

interactions, where M is the number of modes considered and N the number of 

frequency components. Moreover, for each interaction a wave is produced which 

does not "fit" the guide and requires balancing at 4 = 0 by a homogeneous field which 

in general contains a full complement of M modes. The total number of 

3 2 

contributions to the wave field is therefore 9M N . This is truly forbidding, since 
the use of at least a few frequencies and several modes is essential to a realistic 
investigation of a wave propagating in a lined waveguide. For example, two 
frequency components with two modes apiece would require calculation of no less 
than 288 contributions. Nevertheless, we are only prevented from achieving this 
solution by the prodigious time and effort involved. As a research project in and of 
itself this calculation is well within the realm of possibility. However, the problem 
might be better solved using finite element techniques. 

We are now forced to resort to bold approximations if we want to convey 
any information about the effects of the liner’s nonlinearity on attenuation. We 
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consider the nearly-plane wave mode where a L. and c L,, are both much less 
' mn A mn M 

than unity. In this case, the tangents in Eq. (6.39) can be replaced by their 
argumen.s to give 



(6.45) 


where a 



is the fill ratio. 


It is now possible to solve explicitly for a, b, c, and d, where the lack of a subscript 
denotes the nearly-plane wave mode. In essence we have assumed the duct to be 
acoustically narrow, and we must therefore keep in mind at all times the restriction 
to long wavelengths which we have imposed. If used for wider ducts, the 
calculations give attenuations which are much too large, since we have effectively 
neglected beaming. The expressions for the wave numbers are: 


a 


2 


n^/y (n^/y - q^) 

2 2 
n /y + oq^ 


(6.46) 


and 



n^/y q^ (1 +a) 
_ - 

n /y + oq 
^n 


2 a (n^/Y - q^) 
c = - 


2 /„ 2 

n /y + aq^ 


(6.47) 


(6.48) 


We are particularly interested in the value of -Im(b), because it represents the 
attenuation of the wave as it travels down the duct: 
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Figure 6-9 is a graph of the attenuation in dB/m which is plotted against the 
common logarithm of the normalized liner impedance. The graphs are parameter- 
ized against the fill ratio a . There is clearly an optimum liner impedance for a 
given fill ratio which yields the maximum absorption in the duct. 

The magnitude of the characteristic impedance of the material has been 
shown to increase for an intense wave. We can get an idea of the resulting effect on 
attenuation by considering the pure tone expression for the characteristic impe- 
dance of the material, Eq. (3.41a), where K is based on R2» An alternate expression 
of this is Eq. (B.lla) where the quadratic resistivity model is used. An approxima- 
tion of the attenuation rate in the duct at a given point can be found by using the 
particle velocity magnitude in either of the above equations and reading the 
attenuation rate corresponding to this impedance value in Fig. 6-9. 

It can be seen that increasing the amplitude of the wave does not always 
increase the absorption in the duct. In cases where the small-signal impedance is 
below the optimum, a high intensity wave generally is attenuated more rapidly than 
its linear counterpart. On the other hand, if the small-signal impedance is already 
toe* large, an increase in the amplitude of the wave simply brings a further 
deterioration in duct performance. Although this method of determining the 
attenuation of a duct is not highly accurate, it gives a useful mnemonic for the 
de:»igner. 

A duct should therefore be designed so that the small-signal impedance is 
below the optimum for a given fill ratio (which is usually determined by other 
physical constraints). The degree to which the impedance lies below the optimal 




value is determined by the initial intensity of the wave and the length of the duct, 
so that the absorption hovers in the vicinity of the maximum. An algorithm is 
therefore needed which integrates the total attenuation over a given length of lined 
ducting (since the absorptive power of the duct changes in space because of the 
decaying wave). We suggest a procf'dure which calculates the wave amplitude and 
corresponding liner impedance for successive small steps through the waveguide, and 
integrates the corresponding infinitesimal attenuations. The optimal valve would 
then have to be found on a trial and error basis, since no closed form optimization 
formula has been obtained thus far. 


E. Conclusions 


Calculation of the dispersion and nonlinear effects caused by the material 
are complicated in vector notation, but addition of the inherent dispersion of a 
waveguide makes a perturbation solution extremely tedious. We have, however, 
managed to salvage an understanding of the gross effects of the nonlinearity on the 
absorptive power of the liner. Further study is recommended into solution oi the 
modal wave number equations, and computational perturbation or finite-difference 
solution of the nonlinear problem. 


141 


APPENDIX A 


THERMAL EFFECTS ON SOUND SPEED 


This Appendix contains an investigation of the frequency dependent 
effects on sound sp>eed caused by heat conduction between the air and fibers. The 
bulk of this discussion parallels that given by Kuntz^ in an appendix to his Ph. D. 
dissertation. Two assumptions are made: the wavelength is much larger than the 
fiber diameter, and the fluid flow around the cylinder is laminar. The wavelength 
assumption applies well beyond the range of audio frequencies and means that tne 
pressure forcing function is uniform throughout the region of interest. The laminar 
flow assumption breaks down between 125 dB for fl=.98 and 135 dB for R=.94 batted 
Kevlar' 29. Above this level the flow becomes turbulent and the primary mode of 
heat transfer is convection rather than conduction. In a turbulent flow the 
prevailing thermodynamic conditions are expected to be more isothermal than 
predicted for a laminar flow model. 

The flow of heat to and from an individual cylindrical fiber is governed by 
the radial heat conduction equation 



1 9l 1 

r 8r " a at ~ " k 


exp(jcjt) 


(A.l) 


diffusivities. The subscript f refers to quantities pertaining to the fiber, while the 
lack of a subscript indicates a quantity for air. Since the fiber is considered 
incompressible and the wavelength large the total pressure is independent of range 
and can be expressed as 2 =P exp( ja>t)+p^. In this case the radial heat equation is 
solved by 
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T = ( A I (rt) + B K (rt) - p/p c ) exp(jwt) , (A.2) 

— o o o p o 

vk-here K and I are modified Bessel functions of the first and second kinds, 
o o 

respectively, and t =yjw/of . The specific heat Cp should not be confused with the 
phase speed presented in Chapter 3. Since the thermal properties of the air and 
fiber are different we define two temperature functions, T=Texp(ja>t)+T^ for the air 
and Tj=Tjexp(jart)+T^ for the fiber. The subscript f will henceforth indicate a 
quantity appropriate for the fiber, while the lack of a subscript indicates a quantity 
fir air. The temperatuic and heat flux are balanced at the wall of the fiber, r=a, 
and the following temperature distributions are found. 


T = 


o p 


- 


K> t) 
o 

D 


exp(jwt)|+ T 


(A. 3) 


r > a 


where D = K (at) - 
o 







9 


and 



r,- / exp(jwt) + f 


(A.4) 


r<a 


The value of decays at large ranges, and 1^ decays as the range is reduced to 
zero. Temperature fluctuations are therefore largest far away from the fiber and 
are usually near zero at the center of the fiber. Th*.* region in between these ranges 
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contains two thermal boundary layers inside which the transition is accomplished 
and the thermal losses ti.ke place. The boundary layer outside the fiber has a 
thickness 6=y/2aJ(j and the layer inside the fiber, The internal boundary 

layer is much thinner than the external one and can usually be neglected. It turns 
out that the thickness of the external boundary layer plays an important role in 
determining the sound speed because it influences the overall thermodynamic state 
of the fluid in the pores. 

We now evaluate the average temperature fluctuation throughout a 
representative region of a fluid which is filled with a multitude of individual fibers. 
For this analysis we need only assume that the fibers are cylindrical and parallel. A 
spatial average of the value of T is taken over differential concentric circles 
between the fiber wall and some as yet unspecified integration radius r^, 


<T> = -^(1 -<F>) , 

o p 


Ko(rf) 

where F = — — and 


<F> = 


1 


2 2v 

T(ro-a ) 


r 


27rrK (rt)dr 
o 

D 


The result of the integral is 


<F> = 


2 /^KiV)-aKj(at) 


?(r2-a2) 


D 


(A.5) 


(A.6) 


(A.7) 


The average value of F is an indication of the overall thermodynamic state of the 
fluid in the pores. When <F> = 1, conditions are isothermal and the temperature 


fluctuations are very small. When <F> is near zero, adiabatic conditions prevail and 
the effect of the fiber is negligible. 

The integration distance is chosen so that the total integrated volume 
around all of the fibers is equivalent to the volume of fluid in the pores. It was 
found in Chapter 2 that the volume and area porosities are equivalent; the 
expression for the porosity can be written as 

fi=l-Cjra^ , (A.8) 

2 

where C is the number of fibers per unit area and Tra the area of a cylindrical fiber. 

A circle of radius r^ is centered on each fiber, so the percentage of the material 

2 

volume falling within the integration rav^ ■J. .. Cirr^ . Upon eliminoting C between 
the two equations, one finds that 



(A.9) 


Although not all of the fluid parcels are counted, a number of parcels are counted 
more than once to make up the difference. The resulting volume is presumed to be 
a good representative of the fluid. 

The value of <F> is now used to determine the average sound speed in the 
pores. The ideal gas law is written in modified form using spatial average values as 

<£> = <^>R<T> , (A. 10) 


where R is the difference of specific heats Cp-c^. 


Is equation can be written as 


P + 


<p>exp(jwt) = +<p>j exp(jwt)j R |t^ + <T>exp(jwt)j , (A. 11) 
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where the total density has been expanded as were the total pressure ana 
temperature. The "ambient" terms cancel out (p^= and we neglect the 

second order term. The result is 

<p> =R(<p>T^+ P^<T>) . (A. 12) 

Since the expression for <T> is already known, we substitute it .nto Eq. (A.12) to 
find 


RTo 

R7cp(f:<F>r • 

We define the average sound speed <c> as>^p>/d<p>. Since the term in brackets 
in Eq. (A.12) is simply a constant, the differentiation is simple and <c> is clearly 

c 

<c> = ^ , (A.U) 

v/1 + (Y-1)<F> 

where V = 1.4 (for air) is the ratio of specific heats c /c . This is nearly the result 

p V 

obtained by Kuntz;^ the only difference is the choice of integration radius. A 
general discussion of the results of this equation has already been given in 
Chapter 2, B.4. 

In this analysis it is assumed that the external boundary layers of adjacent 
fibers do not interact. At very low frequencies, however, the boundary layers 
become very large. In this case the temperature solution, on which the calculation 
for <c> is based, is no longer valid. We define a criterion frequency, below which 
the calculated value from Eq. (A.14) cannot be trusted, as the frequency at which 
the boundary layers of the fibers are predicted to intersect the centers of adjacent 


fibers. The fiber radius a is usually known for a given material; for Kevlar it is 
6 • 10”^ m. The expression for 6 has been given previously and is ^2afu where 
Q'=2.2 • 10 m /sec is the thermal diffusivity of air under standard conditions. To 
find the inter-fiber spacing, we return to hq. (A. 8), where the average number of 
fibers per unit area,^yc, is defined. The inter-fiber spacing is simply 1/JC, the 
"distance per fiber", 



(A.15) 


Our criterion can be written as s=a+6. We first solve for the critical 6 at which the 
criterion is reached. 



We then solve for the angular frequency 


(A.16) 


2a 



and the physical frequency f^^ below which the solution is inaccurate. 




(A.17) 


The actual sound speed in the pores is expected to be closer to isothermal than 
would be calculated from Eq. (A.13) for frequencies below f^^. For the three 
porosities of Kevlar used in our experiments, the frequencies are 


n=.98 1462 Hz , 


12 =.96 3147 Hz , and 

«=.94 5002 Hz . 

Figure ^-1 is included for the reader's convenience and is merely a reprint of 
Fig, 2-2. It shows that the sound speed tends towards the isothermal value in the 
middle of the audio band. The actual sound speed value is, however, closer to the 
isothermal value than predicted by Eq. (A. 14) for frequencies below f^^. Moreover, 
at high intensity the sound speed is closer to isothermal fo all frequencies because 
of turbulent mixing in the fluid. Since we are most interested in measurements 
taken below f^^^ and above 135 dB, we are justified in assuming ttie sound speed to 
have the isothermal value. 
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Figure A-1 

Frequency-dependent sound speed 



Appendix B 

RELATIONS FOR QUADRATIC RESISTIVITY 


2 

The quadratic resistivity function r=rj+r^u was not used in the bulk of 
this investigation because the rj+r 2 usgn(u) model was observed to fit the static 
flow resistivity data better. The two models, however, closely approximate one 
another for most of the flow velocities used in the static flow resistivity 
measurements. This is to be expected, since both curves are intended to mimic the 
same set of data. Figure B-1 shows corrected measured resistivity data and the 
fitted curves for the two resistivity models. This figure is the same as Fig. 2-1 and 
is presented here for the convenience of the reader. The values of the two 
resistivity functions are close over the measured range. However, the two curves 
rapidly diverge above this range. Table B-1 gives the fitted values of tj and r^ and 
the ratio of the mean-squared errors for the seven samples tested in this study. 


Table 5-1 




'■l 

'■3 


Sample 


(rayl/m) 

2 3 

(raylsec /m ) 

MSE ratio 

1 

0.94 

51007 

4260 

1.04 

2 

0.98 

11214 

1302 

1.86 

3 

0.96 

30700 

4195 

1.18 

4 

0.94 

53031 

8839 

1.42 

5 

0.96 

29497 

4291 

1.38 

6 

0.94 

54394 

6995 

1.40 

7 

0.96 

29727 

3393 

1.70 
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The quadratic resistivity model has one definite advantage; the resulting 
cubic nonlinear term in the wave equation lends itself more readily to analytical 
study than the corresponding u sgn(u) term. We therefore present the following 
results in preparation for a planned perturbation solution for a high intensity sound 
field in a lined duct. 

The Helmholtz relations based on the quadratic resistivity model have the 
cubic not .linearity form 

-nflR.. 

V + jq V = ^5 ^ > V V V 

n '^n n 2q^ / ^ / a P q n-p-q , (B.l) 

P q 

2 

where ^ 3 =*' 3 Cj/ • If the harmonics of the signal are assumed to have a 
negligible effect on the fundamental, we assume that only Vj and Vj* are important 
terms in the convolution. Equation (B.l) becomes 


which reduces to 


jq,V| =- 


3BR, 


(B.2.b) 


This equation can be written in A, 0 notation (see Chapter 3) as 


A'j + j(qj 


P’vi)Ai = 





(B.3) 
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Once again the real ' id imaginary parts separate to give 


A'l . q,|A, = -Kjqi^Aj 
and 


(B.4) 


*^Ir " ®vl ~ ”*^3^1i'^I » 

where 


(B.5) 




J'/T7 


R 


(B.6) 


The solutions of this equation are similar in form to the ones generated for the other 
model .n Chapter 3. 


A^(X) = 


Aj(0) exp(-q^.X) 

\/l +F 3 (i - exp(-2qj.y)j 


where 


K,Af(0) 

r , =- LJ - 


‘^vl'hr'XjqijAj) 


-1 


Z, =q, -jKjAjq* 




-i/.!nlSvN 

qir'XjqiiAj 


2 / ’ 


Zll = (hil^d +k2aJ) + 2RjK3aJ) 


1/2 


(B.7) 


(B.8) 

(B.9) 


(B.IO) 


(B.lla) 


% 
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which can be approximated by 


|Zj| a«|cj|(l + K3AJ) , (B.llb) 


where 


There is very little quantitative difference between the solutions of these 
oqjressions and those derived in Chapter 3 if the particle velocity amplitude lies in 
a range where the resistivity functions have similar values. As an example, the 
velocity at which the two curves cross in Fig. B-1 is 1.45 m/sec. Above this 
velocity the value of the resistivity function for the quadratic model exceeds that 
for the sgn(u) model, and the resistivity curves rapidly diverge for high velocities. 
For a 1 kHz tone in R=0.96 batted Kevlar*^ 29, a 1.45 m/sec particle velocity 
corresponds to an SPL of 153.9 dB. Solutions of these equations for signals weaker 
than 153.9 dB (in this particular case) are therefore expected to be nearly equivalent 
to those in Chapter 3. 

The numerical solution of Eq. (B.l) involves a mere two line code change 
in the existing computer program (see Appendix C). The successive approximation 
for the harmonic impedances is less stable in the quadratic resistivity case because 
of the higher power of u involved. The successive approximation algorithm may be 
unable to converge for high intensity signals, in which case the propagation curves 
are never generated. It turns out, however, that the algorithm is stable for virtually 


qir + KjqjjAj - 


2qiiQiC3 exp(-2qjjX) 

Ql2 + (1 - C 3 exp(-2qjjX)) 


t1 


(B.l 2) 




0 -Qpr3 

I + r. 
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all boundary conditions that do not exceed the aforementioned amplitude limit. 

Figure B-2 contains a comparison of predicted ctnd measured harmonic amplitudes. 

The data plotted here are the same as contained in F;g. 4- 12. Since the solutions 

2 

for the two models are nearly equivalent we have no qualms about using r^u in 
place of r 2 U sgn(u) within the stated restrictions on amplitude. 


Theory data 



(aP) IdS 


Distance (cm) 


00 



Appendix C 

COMPUTER PROGRAM 

Computer program PRPKRV2 is a FORTRAN 5 program which solves 
Eq. (3.29) as an initial value problem. A resistivity function of the form 
r=r j+r 2 usgn(u) is assumed. The heart of the program is the IMSL routine DGEAR, 
which solves systems of differential equations using a 12th order Adams predictor- 
corrector method, also known as Gear's method. The routine is adaptable to a wide 
variety of problems and is particularly designed to solve stiff equations with 
unstable solutions. Fast Fourier transform and plotting routines which were written 
at Applied Research Laboratories, The University of Texas at Aust;.:, were used 
extensively in this program. 

The calculation retains up to 32 harmonics and accepts inputs of SPL and 
phase for as many. The user must also supply the values of rj and r 2 i the 
fundamental frequency, eind the porosity Q. Hardcopy output is provided which gives 
the dimensionless coefficients Rj and R 2 » the initial dimensionless impedance 
magnitudes for the first six harmonics, and the predicted SPL of the first five 
harmonics out to 20 cm beyond the effective origin in ^ cm steps. A plot of the 
SPL predictions is also supplied, and the user has the option of overlaying measured 
data for comparison with predictions. 

2 

The program PRPKRV3, which uses the resistivity relation r=rj+r^u , 
differs from PRPKRV2 in only two places. These are marked in the code and the 
replacement lines are provided. The iterative approximation for the initial 
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impedances has been observed to be less stable in PRPKRV3 than in its predecessor 
and, as a result, the algorithm may never converge and the solution may never be 
executed. The limits ot the algorithm stability are discussed briefly in Appendix B. 

The computer code contains a running commentary on its function which 
will, it is hoped, help the user understand the operation of this program. 



PROGRAM PRPKRV2( INPUT, OUTPUT, KRUD, PLOT, PROPDAT 
, TAPE4-PRCPDAT, TAPE5=OUTPUT, TAPE6«=KRUD) 
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10 CONTINUE 
RETURN 
END 
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